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Fringe Pattern of an Oscillating Fabry-Perot 
Interferometer 


K. D. Mielenz 


(December 17, 1963 


The dependence of the fringe pattern of an oscillating Fabry-Perot interferometer upon 
the vibration amplitude is discussed experimentally and theoretically. If the fringes of the 
stationary interferometer are wide, the fringe pattern tends to disappear at certain values 
of the amplitude. If the stationary fringes are narrow, a splitting of fringes occurs at certain 
amplitudes. In both cases, the stationary pattern reappears, with reduced contrast, at 
intermediate amplitudes. 


The practicability of using these effects for the measurement of vibration amplitudes is 
discussed. 


1. Introduction 


In a recent paper [1],’ a photoelectric spectrum analyzer was described that consisted of 

a Fabry-Perot etalon with a piezoelectric ceramic tube incorporated into the spacer. If an 

a-c voltage of circular frequency w is applied to the ceramic, the length of the etalon, d, varies 
with time, f, according to 

d=d,- d cos wt. (1) 


The static length of the etalon, d), was 10 cm in one experiment, 8 mm in another. The 
vibration amplitude, d;, is of the order of a few wavelengths. It is proportional to the peak 
voltage applied, V,; the ratio was, in both cases, d,/V;=1.5 mu/v. A 60-cycle a-c voltage 
was used. For both experiments, the green line of Hg 198 was chosen to illuminate the 
interferometer. 

The visual observation of the Haidinger fringe patterns of these oscillating interferometers 
reveals a peculiar behavior: 

As the vibration amplitude increases, the fringe pattern of the 10-cm etalon almost dis- 
appears at certain discrete values of the amplitude and reappears at intermediate values; 
figure 1 


With the S-mm etalon, on the other hand, the pattern also reappears at certain amplitudes 
but, instead of vanishing at intermediate amplitudes, it appears as a double pattern; i.e. 
each fringe split into two, figure 2. 


, with 


The contrast of the pattern in both cases gradually decreases as the amplitude increases. 

Fringe disappearance, as in the first case, is well known from oscillating two-beam inter- 
ferometers. Described by Osterberg more than 30 vears ago [2], it is now the basis of an estab- 
lished technique of measuring small vibration amplitudes [3]. 

As of this writing, the obviously more complex behavior of the oscillating multiple-beam 
Fabry-Perot interferometer, however, has not yet been explained. 


2. Theory 


Consider an etalon of length d, and let 7 and FR be the transmission and reflection coeffi- 
cients of the etalon plates. Its transmission is 


A(v, d, 1) =T?(1+ R?—2R cos ¢)=", 
with 
By, 
(4r/c)d cos 7, 
speed of light, 
frequency, 
angle of observation. 


——— 


Figures in brackets indicate the literature references at the end of this paper 
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f i 0, 6 ( 54, and 201 v. 
F 1 F patter of oscillating 10-cm etalon, for \ ) 68, 110, 15 
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Figure 2. Fringe pattern of oscillating 8-mm etalon, for Vi=0, 46, 96, 148, and 178 v 
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Div) Ce~= &—)* Ce a/B)*(" : (5) 
be the Doppler distribution of the illuminating light, ») being its principal frequency. Here, 


C, a=constants, (6) 
Q Bro. 


Monochromatic light, of frequency », is then observed in a given direction 7 with the intensity 
Td, 4 D(v) Av, d, 7). 


The total intensity observed in that direction is 


. 


I(d, ) ) Td. i)dv. 


which, upon substitution of (1), vields the instantaneous intensity distribution in the fringe 
pattern of the oscillating etalon. The actually observed time average, then, is 


I (dy, d,, 7) (1/27) is I(do+d, cos wt,i)d(wt). 


It is probably impossible to solve this integral in a closed form. A solution can be obtained, 
however, by expanding the integrand into a Fourier series as follows. 

In the usual derivation of the Airy formula (2), the complex amplitude of the wave trans- 
mitted by the etalon is obtained as [4] 


P=uT >; 8 
The intensity is 


A=PP*=T? DY RiRte~ ee 


_ T? ( y R*+35 (¢ ing _! ¢ ing) * 2 prt ) 


n=0 n=1 m=0 


— T?/(1—R?) ( 1+2 >> R" cos ng ) 


n=1 


Thus, 


I(d, :)=CT’ i—r)| [ e~(aiB"e-¢ “d(o/B) t 25° R" e~ (a B(e-o “ cos ne ae/a) | 


n=1 0 
» 9 


Since e~‘* **" is essentially zero for — © <r<—@, we have 


r/p 


¥-*” eosne de=cos nd ( e~ (2®)"" cos na dz 
. — & 


—sin nod [- e-(@ B= sin nz dr=cos no(ymB/a)e~ ‘2 "—0. 
Therefore, 
I(d,ij=K (1 +2 > Rre~("820)* eos np )> 
n=1 
with 


K=y)x/a-CT? (1—f?). 
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For the oscillating etalon, 


o=(42v/c) (do +d, cos wt) cos 7 


do +d; COS wt, (12) 
and - 


B/2a= (2m/ac)(d)+d, cos wt) cos 7. (13) 


Neglecting the time as well as the angular dependence of the exponential function, we obtain 


I (do, dy, 1) =1 (bo, 6) =K 1 t 2>°R* 2nxdy/a)” eos (ndo+nd, COS wt | (14) 
n=1 

where ¢ determines the static pattern, and ¢, describes the vibration amplitude. 

Fourier expansion of the integrand in (8). 


Hence, 


This is the 


2° 


I(», $1) Kn) | U (wt) +235 Rre~ 2nxto/a0)* 
JV n=1 


COS (Ndy-+ Nd; COS wt) dot) } 
With 


r °%2r 
COS( Nd) + nd, COS wt) d(wt)=cos ndp | cos(nd, cos wt) d(wt) 
PT All 


Ps 3 


—sin nd, | sin(nd@, cos wt) d(wt) =2r cosndyd (nd) +0, 


J 0 


where J, is the Bessel function of zero order, we finally obtain 


I(¢, 91) > K[ 1 r2 y+ tg —2nxdy/ac)*eosnd, J,(nd, | (15) 


This is the desired time-average of the intensity distribution in the fringe pattern of the 
oscillating etalon. 


3. Discussion 
3.1. Calculation of 27i/,/ac 


For the experiments described in this paper, the green line of Hg 198 was used. 


Its 
half-width is, approximately, 2(v—pp),/2=8-10* see~!. Thus, from (5), 


a’16-10'® see 


and 
2rd, /ac=d, - 
with d) in em. 


3.2. The Long Etalon: Fringe Disappearance 


For the experiment of figure 1, the data of the etalon were d,=10 em and R 


0.8. Thus, 
from (15) and (16), 


n=1 


<I(¢,¢:) > KI +2 S* (0.8)"¢ odd COSN dy Iu(nds) 


This series converges very fast. To obtain <J(@, )) >/<1(0, 0) 


- to better than four decimals, 
it suffices to consider the first two terms of the sum only; n <2. 
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A plot of <J> versus ¢, is shown in figure 3 for ¢=2m7 (bright stationary fringes), 
)-—-(2m+1)mr (dark stationary fringes), and @=(2m+4)r; m=0, +1, +2,.... The 
first two curves show maxima, or minima, at ¢,=0, 3.9, 7.1,... . For intermediate ampli- 
tudes, @,—2.4, 5.5,. . ., the three curves practically cross each other in one point. 

The intensity distribution in the fringe pattern (</> versus ¢») is shown in figure 4 for 
these significant values of @. The pattern practically disappears for ¢,—2.4 and 5.5. It re- 
appears as a “negative” (dark fringes instead of bright ones) for ¢,=3.9, and as a “positive” 
(fringes in their original positions) for ¢,=7.1. As ¢, increases, the contrast decreases. 

The peak voltages V, at which these conditions were actually observed are listed in table 1; 
the corresponding vibration amplitudes d, are proportional to these voltages. As must be 


expected from (12) if the theory is correct, the ratio ¢,/V; is a constant, here equal to 0.035 v™. 




















—- <1(%9,%,) >/<1(0,0) > 


FIGURE 3. versus $, for d)>=2mz, (2m+}%)m, and (2m-+ 
10=cm etalor 























4r 57 
—- 


Figure 4. Theoretical intensity distribution in 


fringe pattern of 10-cm etalon, for ¢,=0, 2.4, 
3.9, 5.5, and 7.1. 


78 





TaBie 1. Significant amplitudes, 1, and peak voltages, Vi, 
for 10-cm etalon 


Pattern p 1 i] Vi 


Dis ippeared 
Negative 
Disappeared _ 
Positive 


3.3. The Short Etalon: Fringe Splitting 


For figure 2, an etalon of 8 mm length and, again, 80 percent reflectivity was used. There- 
fore, from (15) and (16), 


<I(¢,¢1) > K1 +2 >) (0-8)" én? .10~* eosnd, Jo(ne) | 


n=l 


The convergence of this expansion is rather poor. It takes 10 terms of the sum; n<10, to 
obtain <1(, d,) >/<1(0, 0) >to about two decimals. This accuracy, however, is sufficient for 
our purpose. 


The plot of <J> versus @, again for ¢)=2mz, (2m+1)z, and (2m+'s)z, is shown in 
figure 5. Quite contrary to the corresponding graph in figure 3, the three curves cross each 
other in no single point. The fringe pattern can therefore no longer be expected to disappear 
at certain vibration amplitudes. The curves for ¢=2mz and (2m+1)z are no longer mirror 
images of each other; they have maxima at ¢,=0, 6.4,.. ., and ¢,=3.4,. . ., respectively. 


The curve for ¢)>=(2m-+ 4) x, practically a straight horizontal line in figure 3, now has pronounced 
peaks also; at $,=1.7, 4.9, 





- <I(®9,®,) >/<1(0,0) > 











Figure 5. <I >versus ¢; for ¢>=2m7, (2m+4)z, and (2m+1)zx. 
8=mm etalon. 


79 




















>/< (0,0) > 


<1(05,%,) 














2r 3 
~ 


Figure 6. Theoretical intensity distribution in 
fringe pattern of 8-mm etalon, for @,;=0, 1.7, 3.4, 
1.9, and 6.4. 


Figure 6'shows </> versus @», the intensity distribution in the fringe pattern. One sees, 


as in figure 2, (a), the stationary pattern at ¢,=0, (b), the split ‘positive’ pattern at ¢,=1.7, 
(c), the “negative” pattern at ¢,=3.4, (d), the split ‘‘negative’’ pattern at ¢,—4.9, and (e), the 
reappearance of a “positive”? pattern at ¢,=6.4. Again, the overall contrast decreases as 
¢, increases. 

The peak voltages V, at which the photographs in figure 2 were taken are compared to 
these ¢,-values in table 2. Again, the agreement between experiment and theory is satisfactory 
as o,/V, turns out to be the same constant as in table 1. 


TABLE 2. Significant amplitudes, ¢:, and peak voltages, V 
gn, impe } { 
for 8-mm etalon 


4. Conclusions 


The fringe pattern of oscillating Fabry-Perot etalons, as shown in figures 1 and 2, are fully 
explained by the above theory. Fringe disappearance, such as also obtained with an oscillating 
two-beam interferometer, is observed when the fringes in the stationary pattern are broad. 
Fringe splitting is found when they are narrow. The vibration amplitudes significant for 
fringe reappearance, and disappearance or splitting, are different in the two cases. 

The behavior of the pattern, thus, depends entirely upon experimental circumstances 
such as, the length of the etalon, the reflectance of the etalon plates, and the width of the 
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spectrum line used for illumination. Mathematically, these factors determine the rate of 
convergence of the series expansion (15) and, thus, the distinction between the two border cases 

For an oscillating two-beam interferometer, no such dependence upon experimental 
circumstances exists. For the measurement of vibration amplitudes, therefore, an oscillating 
multiple-beam interferometer appears to be less practicable than a two-beam interferometer. 
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Interferometer for Testing Extended Surfaces Such as 
Surface Plates and Precision Ways 


J. B. Saunders and John V. McDermott 
(January 21, 1964) 


A modification of the previously described ‘Interferometer for Large Surfaces” is de- 
scribed that permits the use of smaller prisms for a given sensitivity and a more practical 
arrangement of the optics. This arrangement allows for better structural connection of the 
elements into a rigid and moveable unit. This is obtained by a 90 degree rotation of the 
prism relative to the plane of incidence of the light. 


Since publishing a description of the interferometer 
for large surfaces! these authors have considered 
means of facilitating its construction and appli- 
cations. The modified version permits the use of a 
smaller prism for attaining the same sensitivity. 
The elements are more easily combined (mechani- 


emerge parallel and normal to the prism face. One 
component beam, after emergence from the K6sters 
prism, is reflected normally from a mirror so that it 
returns along its path of incidence to the dividing 
plane. The other component beam is deviated by a 
plane wedge so that it is incident onto the surface to 


be tested. A second mirror receives this beam and 
returns it (after a second reflection from the surface 
being tested) along its previous path to the beam 
divider. The two component beams recombine at 
the beam divider and after a second division a 
component of each beam is received by the observer. 

If the two paths are approximately equalized 
interference fringes of good contrast are observed. 


cally) into a rigid unit. 


Figure 1 is a view of the optical elements. Light 


from a monochromatic source is collimated before 
entering the plane entrance face of a Késters prisin. 
The beam is divided by amplitude division at the 
After total internal 
reflection inside the prism the two component beams 


semireflecting dividing plane. 


MIRROR (MAXIMUM REFLECTIVITY) 


POINT OF INCIDENCE 


KOSTERS PRISM 


EDGE OF SURFACE 
TO BE TESTED 


MIRROR (PARTIAL REFLECTING) 


FIGURE 1. A plan view showing the arrangements of elements 


for a compact surface plate interferometer. 


1 An interferometer for testing large surfaces, by James B, Saunders and Franz 
L. Gross, J. Res. NBS 62 (Apr. 1959) RP2943 
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Figure 2. Interference fringes with light reflected from a 24 


in. black granite surface plate. 


The width and orientation of the fringes may be 
altered at will by adjusting either of the two mirrors. 
Figure 2 show fringes from a 24 in. surface plate of 


black granite. 

Since the reflectivities of all surfaces to be tested 
are well below unity one mirror must have a higher 
reflectivity if the two recombined beams are to be 
approximately equal. 


When testing large surfaces of the type indicated, 
one is usually concerned with measurement of the 
surface along a line. Consequently, the horizontal 
width of the aperture (or field of interference) can be 
relatively small. The height (fig. 1) of the prism 
determines the length of surface that may be tested 
with a given sensitivity (angle of incidence). Thus, 
a prism with a horizontal cross section that is small, 
relative to its height, can be used. This reduces the 
difficulty of constructing a prism having a given 
vertical aperture. 

It will be noted that all elements of this inter- 
ferometer are above the plane of the test surface. 
The mirrors are placed laterally to each other instead 
of one vertically above the other, as in the previously 
described instrument (see footnote 1). These dif- 
ferences allow for simplification of the mechanics of 
the instruments. The sensitivity may be altered by 
replacing the plane wedge with another having a 
different deviation and then readjusting the distance 
between mirrors and prism. 

The sensitivity of this interferometer is the same 
as that described (see footnote 1). One fringe corre- 
sponds to a departure from flatness of \/ (4 cosa), 
where \ is the wavelength and a@ the angle of inci- 
dence. 

(Paper 68C2-153 
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Ferrimagnetic Resonance in Polycrystalline Ferrite and 
Garnet Disks at L-Band Frequencies 


W. E. Case, R. D. Harrington, and L. B. Schmidt 


(January 15, 1964) 


An experimental study made on the measurement of the gyromagnetic ratio and line 
width of four ferrimagnetic materials at L-band using thin disk-shaped samples in a rec- 


tangular TEjo2 mode cavity at 1107 Me/s is described. 


The study emphasized ferrimagnetic 


resonant measurements on a one half inch diameter sample of a given material as a function 


of aspect ratio. 


Experimental results indicate that for three materials the gyromagnetic 


ratio could be accurately calculated from Kittel’s equation for aspect ratios of 50 or greater. 
Because of an anomaly, the calculation was impossible for a fourth material, a nickel alumi- 


num ferrite. 


Generally, very high aspect ratios are required to cause line width measure- 


ments to rapidly approach an asymptotic value. 


1. Introduction 


Considerable effort has recently been directed 
towards the development of ferrite and garnet de- 
vices such as isolators and circulators at the lowe 
microwave frequencies. However, there is a lack of 
engineering information available on the line width 
and gyromagnetic ratio of ferrites at these frequen- 
cies even though these quantities are basic to the 
design of microwave devices. One of the simpler 
methods of measuring line width, A//, and the gyro- 
magnetic ratio, y, at higher microwave frequencies 
such as X-band has been based on the determina- 
tion of the absorption of a small spherical sample in 
a nondegenerate rectangular mode cavity [1, 2].’ 
This technique can be extended to L-band frequen- 
cies although in many cases it may be necessary to 
use thin polycrystalline disk-shaped samples rather 
than spheres in order to avoid magnetic saturation 
problems [3, 4]. However, in the case of disks, it 
is necessary to apply demagnetization corrections in 
calculating y, as was pointed out by Kittel [5]. 
Furthermore, at L-band frequencies, it is necessary 
to use a disk orientation in the cavity which is more 
susceptible to dielectric effects than is usually the 
case at higher frequencies where smaller cavities of 
different modes may be used. It was thus the pur- 
pose of this study to determine the effects of sample 
size and geometry upon the accuracy of L-band fer- 
rimagnetic magnetic resonance measurements. This 
involved an experimental study of line width and 
gyromagnetic ratio as a function of aspect ratio 
(diameter/thickness) of several ferrites. 


2. Experimental Techniques 


The development of a ferrimagnetic resonance 
measurement facility at L-band frequencies is sub- 
ject to two basic requirements. First, the d-c field 


Figures in brackets indicate the literature references at the end of this paper. 


must be perpendicular to the plane of the disk- 
shaped sample and the RF field in order to avoid 
low field losses [4]. Secondly, the cavity must be 
designed so that it can be placed between the pole 
caps of a typical electromagnet. This latter require- 
ment becomes somewhat more severe at L-band due 
to the large size cavity that must be utilized. How- 
ever, both requirements are readily satisfied by plac- 
ing the disk-shaped sample at the center of a TE 
rectangular transmission cavity (fig. 1), as 
pointed out by Sirvetz and Zneimer [3]. 

The measurement procedure utilized for obtaining 
data in this study is similar to an attenuation tech- 
nique which is being issued as an ASTM tentative 
method of test for Ferrimagnetic Resonance Line 
Width and Gyromagnetic Ratio of Nonmetallic 


4102 


was 


DISK 
SAMPLE 


ADJUSTABLE 7 
COUPLING j 
LOoP = TE, 


02 CAVITY 


| RF MAGNETIC LINES 


a 


FicureE 1. Rectangular transmission cavity using TE\o2 


at 1107 Me/s. 


mode 





Magnetic Materials designated ASTM (C524-63T. 
Our setup varies from the ASTM method in that the 
input power to the cavity held constant and 
measurements are made with a precision attenuator 
located in the IF detection circuit following the 
cavity. The frequency is adjusted to give cavity 
resonance for the empty cavity, and the variable 
precision IF attenuator is set to Ae to give a con- 
venient output reference level, as indicated on the 
output detector. The sample is then inserted, the 
d-c magnetic field adjusted for maximum absorption 
in the sample, the frequency adjusted for cavity 
resonance, and the precision attenuator set to detain 
the same output reference signal. The new attenua- 
tor reading is defined as A,, and the d-c field reading 
is defined as H,. The attenuator reading correspond- 
ing to the half-power points on the ferrimagnetic 
resonance loss curve of the sample is calculated from 
the following formula [6]: 


is 


1o—A, ; 
20 +1). 


The attenuator now at the calculated <A, 2 
value. The d-c magnetic field necessary to return 
the output power to its previous level on one side of 
resonance is then read and recorded. The d-e field 
is then similarly adjusted and recorded for the other 
side of The difference between these 
two d-c field measurements is the line width, AH. 

The effective gvromagnetic ratio for the case of a 
sphere is determined from the formula 


A 


A +20 log 2 20 loz (10 


is set 


resonance. 


on f 


afi 


H. (1) 
where /, is the microwave frequency required for 
cavity resonance at maximum sample absorption _— 
IT, is the corresponding d-c field. In the case of ¢ 
disk magnetized perpendicular to its plane, the 
gvromagnetic ratio is given by 


a f 


N 


(2) 
H, »M, 

which is Kittel’s equation where the demagnetizing 
factors N, and N, are equal and N,+N,+N, 
M. is the magnetic moment per unit volume. 

As can be seen in figure 1, the sample is oriented 
with its plane perpendicular to the E-field in the 
cavity. Even though it is at an electric field node, 
the diameter of the sample must be restricted in 
order to avoid errors due dielectric effects. At 
higher frequencies where smaller cavities may be 
used, it is possible to pick a cavity mode in which the 
plane of a disk may be parallel to the electric field 
in the cavity. In these cases, dielectric effects are 
avoided by having a sufficiently thin sample located 
at an electric field node [7]. The necessity of placing 
the plane of the disk perpendicular to the E-field in 
the L-band cavity is not a serious problem due to the 
longer wavelength involved. This was experi- 
mentally verified by comparing H/, and line width 


4r. 
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data on disks of 1 in. diam with corresponding data 
on the same disk reduced to \ in. diam and ground 
to the same aspect ratio. Resonant field measure- 
ments agreed within 0.3 percent for a given sample 
for three different materials. The dielectric effect 
on H, may thus be considered negligible. For a 
fourth material, a distinct diameter dependence 
anomaly was observed, as will be described below. 
Since the dielectric constant of this material is com- 
parable to that of the other three materials, it is 
believed that this anomaly is not due to dielectric 
effects. In the case of line width measurements, 
the variation between different diameter disks may 
be as great as 5 percent for all four materials. This 
small diameter dependence was not considered a 
serious error since the accuracy of the line width 
measurements due to sample variations, etc., was 
also of the order of 5 percent. Furthermore, it will 
be noted later that much larger errors in line width 
may arise unless disks having large diameter to 
thickness ratios are used. It appears that a % in. 
disk offers a good compromise for keeping dielectric 
effects on line widths low while attaining high aspect 
ratios. 

The samples were supported by a rexolite holder 
at the center of the cavity. The dielectric constant 
of the rexolite introduced negligible error in the 
measurement, was determined from similar 
measurements using a polyfoam holder. 

Possible errors due to wall effect were checked by 
measuring the disk at various distances from the 
broad wall of the cavity. The sample is normally 
placed at % of the distance between the broad walls. 
Placing the sample \ and \ of this distance from one 
of the walls resulted in only negligible change in the 
resonant Measurements. 

Many of the materials designed for use at the low 
microwave frequencies have relatively low Curie 
temperatures which may introduce errors in the data 
due to ambient temperature fluctuations. This 
situation is more serious in disk samples since the 
temperature dependent magnetization appears in the 
gvyromagnetic ratio formula given in eq (2). 
These effects were essentially avoided by stabilizing 
the cavity temperature by means of coils carrying 
constant temperature water. 

Having experimentally demonstrated that the 
above effects do not introduce appreciable errors in 
the measurements, the majority of the effort was 
directed towards a study of the dependence of AH 
diameter ) of the disks. 
thickness 
Early measurements indicated that more consistent 
data could be obtained from a given sample that was 
alternately ground down and measured rather than 
comparing data from several different samples of the 
same material. Initially, ultrasonic machining and 
diamond wheel grinding were used to cut a disk 
sample to a given diameter and thickness. <A small 
rotary grinder using aluminum oxide optical finishing 
powder was then used to grind off small increments of 
thickness allowing many measurements on the same 
sample with a given diameter. In general there 


as 


as 


and y on the aspect ratio ( 





appeared to be no dependence of the resonance 
properties on surface finish for three of the materials, 
as was determined by using powders of different grit 
size. However, in the material noted above which 
showed anomalous effects for H,, a surface finish 
dependence was noted. 


3. Results 


The four commercially available materials referred 
to above and for which the line width, AH, and 
resonant field, /7,, were determined as a function of 
the aspect ratio of the disks have the following general 
composition: 

Sample A 

Sample B 


Yttrium iron garnet 
Magnesium, manganese, 
ferrite 
Sample C—Substituted yttrium iron garnet 
Sample D— Nickel aluminum ferrite. 
All measurements were made at 1107 Me/s. 


aluminum 


3.1. Gyromagnetic Ratio 


Figure 2 shows H, versus aspect ratio for the first 
three materials above. It can be seen that in at 
least two cases an asymptotic value of resonant field 
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Figure 2. Resonant field, H,, versus aspect ratio for 


samples of three materials, A, B, and C. 


was not reached even for aspect ratios greater than 
100, which in this case corresponds to disks less than 
0.005 in. in thickness. This of course indicates that 
in such measurements it is not, in general, reasonable 
to assume that one has an infinitely thin disk such 
that demagnetization corrections are not needed. 
However, these curves essentially express the 
dependence of Kittel’s relation as given in (2) on 
demagnetizing factors which are a function of the 
aspect ratio. The relationship between the demag- 
netizing factors of ellipsoidal samples and_ their 
aspect ratios has been given by Stoner [8]. Thus 
(N,—N,) for a given disk may be readily obtained if 
we assume that the demagnetizing factors of an 
ellipsoid and disk of the same aspect ratio are identi- 
cal. The value of magnetization, M,, may be 
determined from an independent measurement using 
a magnetometer. Thus, knowing , and the value 
of N,—N, as well as H, as given on the curve for a 
given aspect ratio, one can calculate the gyromag- 
netic ratio from eq (2). This should be a constant, 
independent of aspect ratio. When this was done 
for the data shown in figure 2, it was found that the 
variation in y as a function of aspect ratio for ratios 
of 50 or higher was less than 1 percent for each of the 
three samples. This was considered as satisfactory 
agreement between Kittel’s equation and _ the 
observed data. 


The validity of the assumption of equivalent 
demagnetizing factors for disks of the same aspect 
ratio was further substantiated by grinding ellip- 
soids from several disks keeping the aspect ratio 
approximately constant. A negligible effect on the 
measured value of H, was observed. This assump- 
tion thus appears to be valid at least for poly- 
crystalline materials for measurements of this type. 
It has been noted in the literature that this situation 
may not necessarily be true for other types of 
measurements such as the study of magnetostatic 
modes in single crystal disks [9]. 

The above discussion implies that Kittel’s equa- 
tion is quite applicable for predicting the gyromag- 
netic ratio of most samples. However, similar 
measurements taken on Material D, which is a 
nickel aluminum ferrite, showed an anomalous be- 
havior as can be seen in figure 3. This material 
varies from the previous ones described in that H, 
decreases for the higher aspect ratios and the data 
is strongly diameter dependent. Recent work in- 
dicates that H, is also dependent upon surface finish. 
The cause of this behavior is not yet understood. 
However, it is apparent that Kittel’s equation as 
given in (2) is not applicable for describing this 
material. It thus appears that care should be taken 
in applying Kittel’s equation to an unknown sample 
without first determining that the material is not 
following an unusual behavior. This implies that it 
may be necessary to measure samples of several 
different aspect ratios when gyromagnetic ratio 
data is desired. 
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3.2. Line Width 


Figure 4 shows the line width of materials A 
through D as a function of the aspect ratio of the 
disk-shaped samples. Here again, it is seen that it 
is difficult to obtain disks sufficiently thin in order 
to obtain an asymptotic value of line width. It 
may thus be necessary to obtain data at several 
different values of aspect ratio and extrapolate the 
results to the higher aspect ratios in order to obtain 
meaningful information on line width. It is readily 
apparent that samples B, C, and D are displaying 
similar phenomena in that they show a relatively 
slow decrease in line width as the aspect ratio in- 
This decrease in line width with increasing 
aspect ratio is probably due to inhomogeneous field 
effects as well as incomplete saturation of the indi- 
vidual crystallites. The inhomogeneous field effect 
might be expected to result in magnetostatic modes 
arising from field inhomogeneities due to spacial 
variations in the demagnetizing field resulting from 
deviation of the disks from true ellipsoids such as 
described by Geschwind and Clogston [10]. The 
problem of incomplete saturation might be expected 
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to cause a line broadening if the relatively low 
internal d-c fields in the disk required for resonance 
at L-band become of the same order of magnitude 
as the anisotropy fields such as discussed by Subl et 


al. [11]. The greatest variation in line width with 
aspect ratio in the data given in figure 4 appears on 
the low field side of the resonance curves. This 
effect might be expected for both line broadening 
mechanisms described above. As such, it does not 
appear possible at the present time to separate the 
eftects due to saturation from the various possible 
inhomogeneous field effects that may be present in 
this data. One further point of interest in regard 
to this data is that sample D does not display an 
anomalous behavior in line width even though there 
was a strong dependence of resonant field on diam- 
eter and surface finish. 

Sample A, which is yttrium iron garnet, shows a 
somewhat different dependence of line width with 
aspect ratio. For this case, the line width appears 
to be relatively constant for the higher aspect ratios. 
As can be seen, it is possible that there is actually a 


slight increase in line width as the aspect ratio 





increases in this material. This apparent increase 
in line width for the higher aspect ratios is not entirely 
understood although here again, the increase is 
largely due to changes in the low field side of the 
resonance curve. It is also of interest to note that 
the resonance curves for this material were in 
general more symmetrical than for the previous 
materials B, C, and D. This is in agreement with 
data obtained on single crystal disks of YIG by 
Smith and Watanabe [12]. The 4 in. diam samples 
of this material proved to be too lossy for measure- 
ments at the lower aspect ratios. However, data 
taken on a smaller * in. diam sample of yttrium iron 
garnet indicated that there is a very large increase 
in line width at aspect ratios of the order of ten or 
less. This is presumably due to the fact that these 
samples are becoming magnetically unsaturated at 
the lower aspect ratios. It should also be noted that 
variations in data on different samples of the same 
diameter cut from the same block were noted for 
this material. This apparent inhomogeneity in such 
materials considerably effects the precision of 
measurements such as described above when different 
samples must be used. 

A final remark concerning the comparison of line 
width measurements between ellipsoids and disks 
of the same aspect ratio for the above four materials 
may be of interest. It has previously been mentioned 
that a negligible variation in resonant field was 
noted when a disk was ground to an ellipsoid of the 
same aspect ratio. However, for the case of samples 
B, C, and D, the ellipsoid had a smaller line width 
than the disk of the same aspect ratio. For sample 
A, the ellipsoid had a greater line width than the 
corresponding disk. This suggests that magneto- 
static mode phenomena affects the line width but 
not the resonant field in these polycrystalline 
materials. 


4. Conclusions 


1. The use of a TE, rectangular cavity is quite 
applicable for measuring the effective line width 
and gyromagnetic ratio of disk-shaped samples at 
L-band frequencies. 


2. In many cases, it may be necessary to obtain 
line width and gyromagnetic ratio as a function of 
aspect ratio and extrapolate the data to obtain the 
asymptotic values for these quantities. However, 


for many materials, the gyromagnetic ratio of a 
finite disk may be calculated from Kittel’s equation 
if one uses the appropriate demagnetizing factors 
and has a value of M available. 

3. In at least one case such as nickel aluminum 
ferrite, it is not possible to fit Kittel’s equation in 
its usual form to the curves giving resonant field 
as a function of aspect ratio. Such an anomaly 
appears to be associated with both a diameter and 
surface finish dependence of H, as a function of the 
aspect ratio. Considerable care should thus be 
given to the interpretation of data of this type on 
such materials at L-band frequencies. 


5. References 

{1] Artman, J. O., and P. E. 

susceptibility tensor in 
No. 9, 1124-1132 (1955). 

Soohoo, R. F., Theory and application of ferrites, pp. 
92-99 (Prentice-Hall, Ine., Englewood Cliffs, New 
Jersey, 1960). 

Sirvetz, M. H., and J. E. Zneimer, Microwave properties 
of polycrystalline rare earth garnets, J. Appl. Phys. 29, 
$31 (1958). 

Lax, B., and K. J. Button, Microwave ferrites and ferri- 
magnetics, pp. 484-485 (McGraw-Hill Company, Inc., 
New York, N.Y., 1962). 

Kittel, C., On the theory of ferromagnetic resonance 
absorption, Phys. Rev. 73, 155-161 (1948). 

Preston, C. C., and W. E., Case, Tables to facilitate the 
determination of the ferrimagnetic resonance line 
width of non-metallic magnetic materials, NBS Tech. 
Note No. 173 (Apr. 1963). 

Von Aulock, W., and J. H. Rowen, Measurement of 
dielectric and magnetic properties of ferromagnetic 
materials at microwave frequencies, Bell System Tech. 
J. XXXVI, No. 2, 427 (1947). 

[8] Stoner, E. C., The demagnetizing factors for ellipsoids, 
Phil. Mag. (7) 36, 803-821 (1945). 

{9] Dillon, J. F., Magnetostatic modes in disks and rods, 

J. Appl. Phys. 31, No. 9, 1605-1614 (Sept. 1960). 

{10] Geschwind, 8., and A. M. Clogston, Narrowing effect of 
dipole forees on inhomogeneously broadened 
Phys. Rev. 108, No. 1, 49-53 (Oct. 1957). 

[11] Suhl, H., L. G. Van Uitert, and J. L. Davis, Ferromag- 
netic resonance in magnesium-manganese aluminum 
ferrite between 160 and 1900 Me, J. Appl. Phys. 26, 
No. 9, 1180 (Sept. 1955). 

{12] Smith, A. W., and A. Watanabe, Ferrimagnetie resonance 
in a single-crystal disk of yttrium iron garnet, J. 
Appl. Phys., Supplement to Volume 32, No. 3, 1555S 
(1961 


Tannenwald, Measurement of 
ferrites, J. Appl. Phys. 26, 


lines, 


(Paper 68C2-154) 








JOURNAL OF RESEARCH of the National Bureau of Standards— 


C. Engineering and Instrumentation 


Vol. 68C, No. 2, April-June 1964 


Changes in the Influence of Atmospheric Humidity 
During Fatigue of an Aluminum Alloy 


John A. Bennett 


(December 


The fatigue strength of 6061-T 
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6 aluminum alloy specimens was found to be about 14 
percent higher in a dry atmosphere than in a moist one. 
after various numbers of cycles 
period during which the presence of water vapor had no effect on the total fatigue life. 


A series of experiments in which the 
indicated that there was an jnitial 
After 


this initial period, cracks developed and propagated much more rapidly in a moist atmosphere 


than in adry one. The change 


in sensitivity to the 


environment is believed to result from 


rupture of the oxide film when the plastic strain becomes locally concentrated. 


I. Introduction 


In recent years there has been a revival of interest 
in the effects of environment on fatigue properties of 
metals. This interest has stemmed from several 
factors including (1) the need for accurate fatigue 
data for aircraft design purposes, so that the effect of 
mild corrosive agents such as water vapor can no 
longer be ignored; (2) the increased emphasis on 
mechanical tests in high vacuum resulting from the 
requirements of the space program; (3) the reali- 
zation that the study of surface reactions can provide 
significant information regarding the mechanism of 
fatigue. 

Results from a number of investigations have 
emphasized the importance of variations in humidity 
on the fatigue properties of aluminum alloys. Liu 
and Corten [1]* found that there was an inverse re- 
lation between the fatigue lives of aluminum alloy 
wires and the moisture content of the air during the 
tests. Although they made no effort to control the 
humidity, the variation in the laboratory over a 
period of months was sufficient to cause a difference 
of more than two to one in the mean fatigue life of 
specimens tested at a given stress amplitude. Coat- 
ing their specimens ‘with petroleum jelly reduced 
the variation, but did not result in as long a life as 
that of the specimens tested in dry air. Broom and 
Nicholson [2] used equipment in which a variety of 
atmospheres, including a high vacuum (2 10~ torr), 
could be maintained around the specimen. They 
summarized the results of their tests on aluminum 
alloys by the conclusion that ‘‘the only constituent 
of normal air that significantly affects the fatigue 
life of age-hardened aluminum alloys is water vapour.” 

The results of the latter investigation provided 
some information on the relative influence of hu- 
midity on crack initiation and crack propagation, 
although there are certain inconsistencies. A few 
experiments in which the environment was changed 


* Figures in brackets indicate the literature references at the end of this paper. 
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during the fatigue test indicated that the total life 
was influenced only by the conditions during the 
latter part of the life, i.e., during crack propaga- 
tion. However, the authors state that these results, 
taken in conjunction with the results of metal- 
lographic examination, “‘confirm that both crack 
initiation and growth are affected by the presence 
of water vapour.”’ 

The observation of gas evolution at the surface of 
specimens during fatigue testing appears to offer 
a new tool for studying the effects of environment, 
particularly on aluminum alloys. As reported by 
Holshouser and Bennett [3], this evolution can be 
seen readily, with specimens of the proper con- 
figuration, if a piece of transparent pressure-sensitive 
tape is applied to the surface of the test section. 
The gas, which is produced only during the latter 
part of the test, forms bubbles in the adhesive and 
lifts the tape from the surface. It was also ob- 
served that the tape had a marked effect on the 
fatigue life, presumably because it restricted the 
access of water vapor to the surface. In a few 
tests the specimens were stressed without tape for 
approximately half of their expected life, then the 
tape was applied; the total life in these tests was 
not significantly different from that obtained when 
the tape had been applied before stressing. 

This result and the observations by Broome and 
Nicholson both suggested that there was a change 
in the nature of the reaction between the metal and 
its environment during the course of the fatigue 
process. In order to verify the occurrence of this 
change and to determine at what point it occurred, 
it appeared desirable to conduct a series of tests in 
which the environment was changed at different 
stages of the test. As the gas evolution phenomenon 
also indicated a change in the nature of surface reac- 
tions, there also appeared to be a need for experi- 
ments to determine if there was a relationship be- 
tween the two changes. 





2. Materials and Testing Technique 


The design of the test specimen was dictated to a 
large extent by the need for a configuration on which 
the evolution of gas could be observed. This re- 
quires that the fatigue cracks initiate in a surface 
to which pressure-sensitive tape can be firmly ap- 
plied, and that the gas be unable to escape from under 
it. Many of the experiments reported in reference 
[3] were conducted with specimens of rectangular 
cross section which were stressed in torsion. These 
met the above requirement, but because of the ex- 
tensive cracking that preceded fracture, it was diffi- 
cult to define failure in a satisfactory manner. Also 
the probability of crack initiation is the same on two 
faces, one of which is difficult to inspect. For these 
reasons it was considered preferable to use specimens 
stressed bending. The tape requirement can be 
met with a bend specimen if the cross section is 
thicker in the middle than it is at the edges so that 
the cracks do not start at the corners of the section. 
The design finally selected is shown in figure la. 

Specimens were cut from a 120° sector of tubing, 
and the reduced section was formed by a 2-in. radius 
cut tangent to the inside surface of the tubing. They 
were stressed in cantilever bending, and it was oer 
that all of the cracks started in the top surface if : 
mean tensile stress equal to '4 of the stress smalls 
was present in this surface. Tape could readily be 
applied to cover the test area, as it was singly curved 
with a radius of * in. 

The specimens were made from commercial 6061 
T6 aluminum alloy tubing. The selection of this 
alloy was based on the results of a previous investi- 
gation [4] in which the fatigue properties of 6061—T6 
specimens were found to show less dispersion than 
either 2024-T4 or 7075—-T6. The aging behavior 
of the material was determined with samples from 
each piece of tubing used for specimen blanks. The 
temperature used (160 °C 


aging was 320 °F 


FIGURE la 


Sian ficant dimensions of the 8 pet 


Maximum hardness was attained in about 30 hr and 
no significant drop occurred in an additional 60 hr. 
Consequently, the specimen blanks were solution 
treated at 970 °F (520 °C) for 30 min, quenched in 
water, and aged 48 hr at 320 °F (160 °C). This 
resulted in an average hardness of about 96 
Rockwell F. 

The specimens were machined and polished after 
heat treatment. Coarse circumferential polishing 
was carried out before the tubes were sectioned. 
After the specimens were machined they were 
polished longitudinally, the final operation being 
hand polishing with No. 600 paper. In order to have 
the abrasive conform to the curvature of the specimen 
adhesive-backed paper was mounted in a tube having 
an inside diameter of 1% in. For the sake of uni- 
formity, all specimens were given a brief final polish 
and inspection on the same day that they were to be 
tested. They were then measured as described below 
and finally they were washed in two successive baths 
of clean benzene. On removal from the second bath 
thev were immediately inserted inthe testing machine. 

The minimum section of this specimen is a seg- 
ment of a circle, so it was necessary to know the 
radius of curvature and thickness in order to calcu- 
late the minimum section modulus. Both of these 
dimensions were measured on each specimen, as it 
was found that the variations in radius were too 
large to be ignored. 

The fatigue testing machine used for these experi- 
ments was of the constant load-amplitude type, 
the load being applied by an eccentric weight rotated 
1800 rpm by asynchronous motor. The alternating 
load was applied to the right end of the specimen as 
shown in figure 1b; any variation in the distance 
from the pivot pin to the minimum section of the 
specimen would, of course, cause a variation in the 
bending moment at that section. ‘onsequently, in 
machining the specimens the holes at this end were 
used to position it for the transverse cut. In mount- 
ing the specimen in the machine, close-fitting screws 
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FIGURE 1b. Photograph showing a specimen mounted in the 


testing machine. 


The alternating load is applied through the connecting rod and pivot assembly 
it the right. Full size 


were used in these holes, so that the moment arm 
was the same for all specimens. In all tests the 
mean load (directed so that the upper surface of the 
specimen was in tension) was one-third of the 
alternating load. 

The entire testing machine was enclosed in a 
cabinet for humidity control. A motion picture 
camera Was mounted above the specimen for taking 
time lapse exposures during tests of taped specimens, 
and a telescope passing through the top of the cab- 
inet permitted viewing the test area at about 20> 
magnification. The humidity control was not auto- 
matic, so the actual values of relative humidity for 
the tests listed as ‘“‘moist environment” varied from 
85 to 95 percent while those for the ‘dry environ- 
ment’’ ranged from less than 5 percent to 9 percent. 
These measurements were made with a commercial 
electrical hygrometer. 

The experiments described in part 4 involved the 
application of transparent pressure-sensitive tape 
to the surface of the specimen. The tape used for 
these experiments was the ordinary cellophane-base 
type; one brand was used for all of the principal 
investigation, but supplementary tests indicated that 
a number of other brands, and even other types of 
tape, would make the gas evolution visible in much 
the same way. 

A few preliminary tests were conducted with taped 
specimens in a low-humidity environment; no bub- 
bling whatever was observed under these conditions. 
Consequently all the “‘taped”’ tests described 
part 4 were run in a moist atmosphere. 

An adjustable contactor was used to stop the 
machine shortly before complete fracture of the 
specimen. Each of the contacts was a 
platinum wire; one of these was mechanically and 
electrically connected to the connecting rod between 
the eccentric weight and the pivot pin. The other 
was on an insulated mount whose position relative 
to the moving contact could be adjusted by a screw 
from outside the cabinet. A dial indicator permitted 
reading the separation of the contacts to the nearest 


piece of 


0.0005 in. After a test was started the contactor 
was raised until it just touched the grounded wire 
at the bottom of the connecting rod stroke, as indi- 
cated by the action of an electronic relay. The 
contactor was then backed off the desired distance 
and the relay set to turn off the machine. It had 
been hoped that this device could be used to indicate 
the development of a small crack by setting the gap 
to a small value. However, the length of crack 
which resulted in cutoff was quite variable and con- 
sequently the technique was not useful in this inves- 
tigation. The contactor did provide a reproducible 
method of stopping the test just short of fracture; it 
was set to a gap of 0.030 in., which resulted in a 
crack through about two-thirds of the section at 
cutoff. This preserved the original configuration of 
the specimen so it could be inspected to see if there 
were any surface features of significance. Also by 
completing the fracture manually after the specimen 
had been removed from the machine, the undamaged 
fracture surface was available for examination. 
This criterion of failure resulted in values of fatigue 
life only slightly smaller than would have been 
obtained at complete fracture. 
It was found that considerable care was required 
in starting the machine to avoid excessive load. A 
variable autotransformer was used to increase the 
alien on the motor very slowly until it pulled into 
synchronism. If this was not done the motor ap- 
parently overshot the synchronous speed, resulting 
in a few cycles of excessive deflection. 
Although a total of three environments were used 
in this investigation (dry, moist, and taped), the 


env ironment was not changed more than once 1M any 


one test. Thus a total of nine types of test were 

conducted, as follows: 

(1) Dry environment throughout 

(2) Moist, then dry 

(3) Dry, then moist 

(4) Moist environment throughout 

(5) Dry, then taped 

(6) Moist, then taped 

(7) Taped, then dry 

(S) Taped, then moist 

(9) Taped throughout. 

The greatest emphasis was placed on tests which 

did not involve the use of tape (types one through 

four), and these will be discussed in the next section. 

2.1. Symbols 

S,—Stress amplitude. 

N,—Number of cycles run in the first environment 
(in tests in which the environment was 
changed). 

Number of cycles to failure. 

In the assumed relations such as those in figure 
5, the number of cycles to the break in the 
curve. 

The estimated number of cycles to form a crack 
0.1 mm long. 

For two-environment _ tests, 
expected on the basis of : 
tion. 

Ny minus N, 


the life to be 
i linear interpola- 





3. Results with Bare Specimens 





Conventional S—N curves for specimens tested in 
both dry and moist environments are shown in 
figure 2. The points represent logarithmic means 
of the individual fatigue lives in all cases, the number 
of replicate tests being indicated on the graph. 
At 26.5 and 38 ksi** the range of plus and minus one 
standard deviation is also shown. The average life 
of specimens tested in the dry environment was 
from two to three times that in high humidity. 
This difference corresponds to a difference in stress 
amplitude of approximately 14 percent throughout 
the range investigated. Although these environments 
were artificially produced, they were only slightly 
more divergent than conditions to be expected in a 
non-airconditioned laboratory during the course of 
a vear. These results emphasize the possibility 
that humidity variations may have contributed to 
the scatter of results in many reported studies of 
fatigue. 

Tests in which the environment was changed dur- 
ing the test were conducted at stress amplitudes of 
26.5 and 38 ksi. The machine was stopped during 
the change, as this required as much as two hours in 
going from moist to dry atmospheres. The results 
of these tests are plotted in figures 3 and 4. Graph 
“A” in each figure shows results for specimens that 
were run for V, cycles in dry atmosphere, then run 
to failure in a moist environment. The abscissas 
are the cycle ratios represented by N,, that is the 
ratio of the number of cycles run in dry atmosphere 
to the mean life in that environment. The “B” 
graphs show similar results for specimens stressed 
first in a moist environment; the abscissas in these 
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graphs represent cycle ratios based on the mean 
life in a high-humidity atmosphere. The end points 
for the two graphs in each figure are the same, as 
they are the number of cycles to failure for specimens 
tested in a single environment. 

As in all fatigue testing, the dispersion of the 
results is troublesome in attempting to determine 
the true nature of the behavior indicated by the 
data. In most cases the logarithmic mean values 
shown in figures 3 and 4 are thought to give reason- 
ably good indications of the true central value; 
however, the point at 33 percent in graph “B” of 
figure 4 is probably unduly influenced by one very 
low result. The mean value of the other two results 
in this group is shown by the cross, and this appears 

; much more consistent with the rest of the data. 
dais: ee name a accel It is apparent from figures 3 and 4 that the fatigue 

ac re life is not a linear function of the cycle ratio in the 
ee ev inticnte the number: of replieete tects, first environment; with the exception of the point 
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single-environment values, and all of the points in 
the “B”’ graphs are above this line. The divergence 
is particularly large in tests run first in the moist 
environment; at a stress amplitude of 38 ksi speci- 
mens stressed for more than 70 percent of their 
expected life in this environment had a mean life 
only a few percent less than that of specimens tested 
throughout in a dry atmosphere. Likewise at 26.5 
ksi the initial stressing in moist environment seemed 
to have little effect on the total life; actually the 
specimen that ran the longest of all those tested at 
this stress amplitude was one that had been stressed 
in moist atmosphere for 30 percent of its expected 
life. 

Even though the dispersion of the data was not 
thought to be excessive, it was large enough to make 
it difficult to fit trend lines directly. Instead it was 
considered more informative to make the simplest 
possible assumptions regarding the behavior of the 
specimens, then determine if the observed data are 
consistent with this hypothetical behavior. As the 
data are clearly not consistent with the linear relation 
in figures 3 and 4, it might be assumed that the 
relation is actually two straight lines, the first one 
horizontal. This would mean that the total life is 
not influenced by the environment during the early 
part of the test. Such hypothetical relationships 
for the tests at 38 ksi are shown in figure 5. If the 
humidity has no effect during the early part of the 
test, then the break in the curves must come at the 
same number of cycles for both environments. This 
number, indicated by Nz on figure 5, was determined 
by trying different values until the best fit for both 
sets of data was found. 

The data were compared with the hypothetical 
relations by plotting the deviations from the linear 
dependence of Ny on N;. The life to be expected 
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FiGuRE 5. Assumed relationship between N; and Ny, for two- 
environmental tests at S4=38 ksi. 


Curve A, dry environment first; curve B, moist environment first. 
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FicurRE 6. Difference between observed fatigue life and that 
expected on the basis of a linear dependence on N,, tests at 
Sa 26.5 ksi. 
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FiGuRE 7. Data similar to figure 6, for tests at S4=38 ksi. 


on the basis of the straight lines in figures 3 and 4 
is designated N,; figures 6 and 7 show the differences 
between the observed mean N; and N, for the same 


values of N;. The solid lines on these figures repre- 
sent the values to be expected on the basis of the 
assumptions of figure 5, i.e., that the total life is 
influenced by the environment only after Nz. Al- 
though many of the points are far from these lines, 
the deviations are both positive and negative, and 
the sums of the deviations in each figure are much 
less than the deviations from a linear relation. Thus 
the assumptions shown in figure 5 appear to be con- 
sistent with the results; in view of the dispersion 
of the data it is probably not worthwhile to see if 
a more complicated hypothesis would provide a 
better fit. Results of the tests with bare specimens 
are summarized in table 1. 





TABLE 1. Results of tests with bare specimens 


Stress amplitude 


A tmosphere 


4. Results with Taped Specimens 


Figure 8 shows enlargements of a few of the time- 
lapse exposures taken of a taped specimen. Prior to 
the application of the tape this specimen had been 
subjected to 900 10° cycles at a stress amplitude 
of 26.5 ksi in a dry environment, but the development 
of bubbles under the tape is typical of all the speci- 
mens photographed. No change is observed during 
the early part of the test, but once the first bubble 
appears the gas continues be evolved until the 
end of the test. In the test of figure 8 the first 
bubble appeared at 1217 10° eveles and failure 
occurred at 1455 10 eye ‘les. 

In the torsion experiments reported in reference 
was found that gas evolution started before 
definite fatigue cracks could be detected either by 
surface examination or by metallographic examina- 
tion of the cross section of the specimen. In the 
current experiments, when tests of taped 
specimens were stopped as soon as gas evolution was 
observed, it was always possible to see small cracks 
under the bubbles. If the tests were monitored 
carefully, the cracks were only au few tenths of 
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FiGuRE 8. Time lapse exposures taken during test of a taped 
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FIGURE 9. S—N curves for specimens tested with tape on the 


surface in moist atmosphere. 


The dashed line is the curve from figure 2 for bare specimens in low humidity. 


millimeter in length and could be observed only on a 
metallographic microscope. The first appearance 
of a bubble on the motion picture film provided a 
convenient and sensitive means of separating the 
fatigue life into two stages. It was estimated that 
the crack length at this point in a test was about 
0.1mm. Consequently, in tests of taped specimens, 
N. was taken as the number of evcles at which ga 
evolution was first observed on the films. In me 
to obtain values of N, for the other two environ- 
mental conditions (bare specimens), a series of 
two-environment tests was conducted in which tape 
was applied or removed. 

Figure 9 shows the S-N curves for specimens 
covered with pressure-sensitive tape throughout the 
test. The open circles represent the logarithmic 
mean life to failure from tests of three or four 
specimens at each stress level. The triangular 
points are median values of .V,, as defined above. 
The dashed line is taken from figure 2 and represents 
the data for bare specimens tested in a dry environ- 
ment. The results with taped specimens differed 
significantly from this curve only at the low stress end. 

Four types of test were conducted in order to 
estimate the values of V, for bare specimens; these 
are listed as types 5 through 8 in part 2 above. In 
types 5 and 6, the bare specimens were stressed for 
N, cycles, then were taped for the remainder of the 
test. If N, was more than N,, gas evolution was 
observed immediately on restarting the testing 
machine. If N, was less than N,, the time lapse 
pictures could be used to determine the actual value 
of N. for the specimen. A sequential procedure was 
used for these tests in an effort to assign values of 
N, close to N,. 

The tests of types 7 and 8 (initially taped) were 
monitored as well as being photographed, and the 
test was stopped as soon as definite bubbles were 
seen. The specimen was removed from the machine 
for cleaning, then replaced and run to failure in 
either a dry or moist atmosphere. 
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The purpose of all types of tests involving the use 
of tape was to provide estimates of N, ‘for bare 
specimens to be compared with the values of Np. 
The tests with initially bare specimens bracketed N, 
but because of the dispersion of the value from 
specimen to specimen this was a very inefficient type 
of determination. The results of tests in which the 
specimens were initially taped gave values of Ne—N, 
(designated N,) for bare specimens. The difference 
between these values and the appropriate mean life 
(Ny) listed in table 1 gave the second estimate of 
N. for each condition. Both of these determinations 
involve some approximations that cannot be rigor- 
ously justified, so the values given in table 2 must 
be considered as only semi-quantitative. 


PaBLeE 2. Results from two-environment tests using tape 


Stress amplitude Atmosphere : Np 


5 ksi Dry 1, 10010 310X103 1, 060 108 
5 Moist 200 120 300 


Dry 37 85 
Moist 2 25 


From tests with initially bare specimens. 
2) Ne=Nr—N; 


In attempting to arrive at a “best value’ of N, 
from the above data, it is not suitable to average the 
values in columns 3 and 5 of table 2. For example, 
the number 22 in column 5 is the difference of two 
much larger numbers, so that its uncertainty is 
greater than that of the other result for the same 
conditions. Consideration of such factors leads to 


the estimates listed in table 3, of the most probable 


values of N,. Data from table 1 


and from figure 9 
are included for comparison. 


TABLE 3. Summary of test results 


Stress amplitude Nr Atmosphere N NP 


210X108 | ~Dry 1, 100108 1, 369X 10 
Moist 250 423 
Taped 5 1, 033 
Dry 
Moist 
Taped 


Although Nz and N, are given to two significant 
figures, it must be remembered that they are based 
on data that show a large dispersion. It is difficult 
to conceive of a mechanism that would result in 
a value of Np larger than N,, as shown for one 
condition, and this is probably not actually the case. 
On the other hand it is believed that the magnitudes 
are approximately correct and that, for example, the 
difference between Nez and N, in the first line of the 
table is certainly significant. Thus the data indicate 
that N, represents a discontinuity in the fatigue 
process which is distinct from the discontinuity 
represented by the development of the first crack. 

These results emphasize the magnitude of the 
influence that water vapor has on the fatigue be- 
havior of this alloy, since the entire effect takes 
place between Np and Ny. At 26.5 ksi this portion 
of the life is five times as long in a dry atmosphere 


as it is in a moist one, and the ratio is four to one 
at S,=38 ksi. Another interesting point brought 
out by these results is that at a stress amplitude of 
38 ksi, the life in low-humidity air of a specimen 
containing a visible fatigue crack is caine than the 
total life of a specimen in a moist environment. 


5. Discussion 


The evolution of gas makes it clear that a chemical 
reaction is taking place during crack propagation 
in this aluminum alloy. It was deduced in reference 
[3] that the gas is hydrogen resulting from a reaction 
between water and the fresh metal surface exposed 
during severe local deformation or crack propaga- 
tion. This conclusion is substantiated by the 
observation that no gas evolution is observed when 
the test is conducted in a low-humidity atmosphere. 

Other observations made during this investiga- 
tion provide some further clarification of the con- 
ditions under which the gas evolution can be observed 
by the development of bubbles under transparent 
pressure-sensitive tape. As mentioned earlier, the 
bubbles on these bend specimens were always 
associated with visible cracks whereas in the earlier 
work on torsion specimens of 5052—-H34 alloy gas 
evolution was observed before cracks could be de- 
tected. This suggests that the development of 
bubbles depends on the rate of gas production; if 
this rate is low the hydrogen can diffuse out through 
the tape before the pressure builds up enough to 
form a bubble. It would be expected that the rate 
would be affected by the area of fresh metal being 
exposed and by the supply of moisture available 
at the metal surface. Apparently the oxide film 
was ruptured over a greater area before crack forma- 
tion in the torsion tests of the cold-worked alloy 
than in the bend tests of the precipitation-hardened 
material. Likewise, when the supply of water is 
reduced by running the test in low humidity, the 
rate of reaction may not be sufficient to result in 
bubble formation. Thus the point at which the 
first bubble is observed probably does not represent 
a discontinuity in the fatigue process but must 
be interpreted as the point at which the rate of gas 
evolution reaches a critical value. 

On the other hand, the results from two-environ- 
ment tests of bare specimens indicate that there is 
a period during the early part of the fatigue process 
when the humidity of the environment has no effect. 
As there is a marked effect after this initial period, 
there appears to be a definite point where the nature 
of the surface is changed so that a reaction with the 
environment can take place. The data of table 3 
indicate that this point is significantly earlier in the 
fatigue life than the start of a detectable crack, at 
least for specimens run in a dry environment. At 
both stress levels the values of Ne and N, for the 
moist environment are nearly equal, while for the 
dry tests there is a large difference in the two values 
at 26.5 ksi, and the difference at 38 ksi is probably 
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significant. These results suggest that the greatest 
effect of the reaction with water may occur between 
Np, and N,; in moist environment a crack forms 
very soon after Np se ican this process takes much 
longer under low-humidity conditions. 

The data of figure 9 show that the tape limits the 
access of water to the surface of the metal to such an 
extent that the fatigue behavior is nearly the same 
as that of bare specimens in a dry environment. 
The evolution of under the tape shows that 
some moisture is reaching the surface, so the results 
of figure 9 indicate that the damaging effect of 
moisture depends on a fairly high rate of reaction. 
This is also shown by the fact that the S-N curve 
for taped specimens coincides with that for dry 
specimens at the high-stress end, but is significantly 
lower at the lower stress amplitudes where more 
time is available for water to diffuse through the 
tape. 

The results of this investigation appear to be 
consistent with the view that the humidity starts to 
affect the fatigue behavior when the deformation 
becomes localized, thus breaking the oxide film. 
During the early part of the fatigue life the deforma- 
tion is apparently sufficiently uniform so that the 
film remains intact and no reaction with the environ- 
ment occurs. Once the film is ruptured, cracks 
develop rapidly if moisture is available, much more 
slowly if the supply is limited. The reason for this 
effect not clear, but it apparently involves the 
oxidation of the aluminum with the liberation of 
hydrogen. 


Fas 


IS 


One possibility that has been suggested 


FIGURE 10. 


[5] is that if the fresh surface formed during slip in 
one direction 1s oxidized, reverse slip cannot occur 
on the same plane, so the damage from repeated 
slip is cumulative. If the freshly formed surface is 
not oxidized, the slip can be reversed on the same 
plane, thus canceling much of the effect of the 
forward slip. 


6. Miscellaneous Observations 


A number of observations that were not related to 
the primary purpose of the investigation are consid- 
ered to be of some general interest. 


6.1. Amount of Shear Plane Cracking 


The fracture surfaces of all specimens were 
examined after testing. As expected, the initial 
plane of the cracks appeared to make an angle of 
approximately 45° with the stress direction, but at 
a later stage the fracture surface was generally 
perpendicular to this direction. All of the fractures 
were irregular, but there seemed to be a definite 
influence of stress amplitude and environment on 
the amount of shear plane cracking. As shown in 
figure 10, the specimens tested at high stress ampli- 
tudes and in dry atmospheres were apt to have a 
large area of shear plane fracture, while in those 
tested at low or high humidity the shear 
areas were found only close to the origin. (The 
section shown in figure 10b did not pass close enough 
to the origin to intersect any such areas.) These 


stress 


Cross sections of fractured specimens intersecting 
the fracture (top) and the upper surface (right). 


50. 


S4=45 ksi, tape on surface throughout test. 


b. S4=22 ksi, moist environment. 
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results are contrary to those reported by Forsyth, 
Stubbington, and Clark [6]; undoubtedly there are 
other factors, such as mean load, that have an 
influence. 


6.2. Color on Fracture Surface 


With many of the specimens, it was observed that 
certain areas of the fracture surface produce bright 
colors when obliquely illuminated with a concen- 
trated light source. By mounting the fracture sur- 
faces in the fractographic specimen mount described 
in reference [7], the colored areas could be observed 
at 50 to 150 on a metallograph. The colors 
appeared brightest when the specimen was. illumi- 
nated directly (rather than through the vertical 
illuminator) with an angle of about 45° between the 
incident light and the axis of the microscope. The 
individual areas showing color were not more than 
a few tenths of a millimeter in extent, but in some 
parts of the fracture a considerable proportion of the 
total area was made up of these colored spots. In 
all cases the greatest amount of color was observed 
in the portion of the fracture where the crack had 
propagated rapidly. No color was observed on the 
static portion of the fracture nor on specimens that 
had been broken in low humidity at relatively low 
levels (conditions resulting in slow crack 
propagation). 

The colors appeared to be spectrum colors rather 
than subtractive colors such as those observed due 
to interference films. Changing the angle between 


stress 


the surface and the direction of viewing resulted in 


a change in the observed color. When the color- 
producing areas were examined at higher magnifica- 
tion with vertical illumination, they were always 
found to contain well-defined striations on the sur- 


Figure 11. Typical electron 


face. This suggested that the observed colors may 
be due to diffraction from these striated areas 

Replicas were made of a fracture surface that 
showed colored areas, and the electron micrographs 
from these replicas (fig. 11) made it evident that 
the striations are sufficiently well defined to act as 
a diffraction grating. The striation spacings meas- 
ured on several micrographs ranged from 0.1 to 
1.3u. As spacings in this range are consistent with 
low-order diffraction of visible light, this is probably 
the correct explanation for the observed color. 

There is considerable current interest in the topog- 
raphy of fatigue fracture surfaces, as evidenced, for 
example, by the excellent investigations of Forsyth 
and others at the Royal Aircraft Establishment. [8]. 
It appears that the diffraction grating effect might 
provide a method of measuring striation spacings 
that could be used for areas that are not accessible 
for direct observation at high magnifications. 


6.3. Usefulness of the Gas Evolution Phenomenon 
in Fatigue Investigations 


The time lapse motion pictures of taped specimens 
that were taken during this investigation have 
suggested a number of ways in which the phenomenon 
of gas evolution might be of use in fatigue studies. 
The bubble formation under the tape provides an 
extremely convenient and sensitive crack detection 
device. As such it has the great advantage that it 
is easily photographed; this permits both accurate 
location of the crack and accurate determination of 
the number of cycles at which the first bubble 
formed. It appears that these capabilities should 
be of value, for example, in statistical studies where 
it is desirable to know the actual stress at the point 


micrographs of a _ fractured 


surface on which colored areas were observed. 
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of crack initiation, or in the many cases when it is 
helpful to separate the fatigue life into periods of 
crack initiation and crack propagation. 


7. Conclusions 


Results of reversed bending fatigue tests on 
specimens of 6061-T6 aluminum alloy appear to 
justify the following conclusions: 


1. The fatigue life in an environment of about 5 
percent relative humidity was two to three times as 
long as that in high humidity. This difference cor- 
responds to a difference in stress amplitude of about 
14 percent throughout the range investigated. 

2. The effect of humidity on the number of cycles 
to failure was much more pronounced during the 
latter part of the test than during the initial portion. 

3. The data are consistent with the view that there 
Was an initial period during which the humidity had 
no effect on the fatigue behavior of this alloy. 

1. The portion of the test during which the 

humidity had no effect was significantly shorter than 
the number of cycles required to form a detectable 
erack in a low humidity environment. In a moist 
atmosphere cracks apparently developed very soon 
after this initial period. 
5. At a stress amplitude of 38 ksi, the mean life of 
specimens containing a small fatigue crack, tested 
ina dry environment, was greater than the total life 
of specimens tested in a moist atmosphere. 

6. Gas evolution under pressure-sensitive tape 
was observed only after visible fatigue cracks were 
present 


7. The development of bubbles under transparent 
tupe provided a convenient and sensitive method of 
detecting small cracks. 

S. No gas evolution was observed when the test 
was conducted in a low-humidity environment. 

%. The extent of shear plane cracking increased 


with increasing stress amplitude and with decreasing 
humidity. 

10. The results seem to be consistent with the 
view that during the early part of the fatigue process 
the deformation is fairly uniform, so that the local 
strain is nowhere sufficient to break the oxide film. 
As stressing proceeds, the deformation becomes 
localized, breaking the film so that the fresh metal 
surface is exposed to the environment. 

11. The striations on some of the fracture surfaces 
were sufficiently pronounced to produce bright colors 
by diffraction. 

There does not appear to be any reason to believe 
that the results of this investigation are unique to 
the alloy investigated; it is probable that they are 
typical of age-hardened aluminum alloys in general. 
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Gas-Lubricated Bearing 
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Iixternally pressurized gas-lubricated bearings have been used successfully to support 


the shaft of a miniature, high-speed expansion turbine. 


cavities to suppress whirl at high speeds. 
detail. 


1. Introduction 


In the field of gas liquefaction and low temperature 
refrigeration, small high-speed expansion turbines 
are often employed as the source of refrigeration. In 
the smaller sizes of machines however, the life and 
speed of conventional bearings are limited, and con- 
sequently a lower limit is set to the practical size of 
such machines. This limitation may be alleviated 
by the use of gas-lubricated bearings to support the 
shaft. The bearings may be lubricated with process 
gas, thus eliminating the risk of contamination with 
oil. Since there is no solid friction, there should be 
no wear and consequently the life of the bearings 
should be of indefinite duration. 

In the operation of a gas-lubricated bearing, the 
surface of the shaft must never touch the surface of 
the bearing, since at the high-surface speeds usually 
employed, such contact may result in seizure and 
the consequent destruction of the bearing, Under 
normal operating conditions the shaft should be 
supported by forces due to pressure in the gas film 
between the shaft and the bearing and these forces 
must be such that the shaft is maintained in a 
position of stable equilibrium. 

In certain types of gas-lubricated bearings, a form 
of instability known as whirl is sometimes observed, 
particularly at high speeds of operation. If the 
speed of the shaft is increased above a critical value, 
the center of the shaft orbits about the center of the 
bearing at a speed approximately one-half the speed 
of the shaft—commonly known as half-speed whirl. 

Whirl may be stable or unstable: Stable whirl is 
characterized by the center of the shaft moving in 
a closed orbit about its former position of equilibrium. 
This is illustrated in figure 1. Unstable whirl (fig. 
2) is characterized by the center of the shaft moving 
along a spiral orbit of increasing amplitude. The 
extent of the increase in amplitude is, of course, 
limited to the extent of the radial clearance between 
the shaft and the bearing. 


The bearings incorporate stabilizing 


The theory of the bearing design is presented in 


The bearing described below has been specially 
designed to ensure stability even at the highest 
speeds of operation with helium gas as the pressur- 
izing medium. The theory of this bearing was pre- 
sented at the First International Symposium on Gas 
Lubricated Bearings [Sixsmith, 1959b]. In the pres- 
ent paper, the theory has been extended and applied 
to the design of the bearings of a helium-expansion 
turbine [Birmingham, Sixsmith, Wilson, 1962]. The 


<_ 


FiGurRE 1. Whirl: constant amplitude. 





Figure 2. Whirl: increasing amplitude. 





theoretical approach used in this particular design 
and construction is outlined in the following detailed 
analysis. 


2. Shaft Whirl 


In the diagram shown in figure 2, it is assumed 
that the shaft is vertical so that the only forces 
acting on the shaft are those due to pressure in the 
surrounding gas. At a given instant of time, the 
centers of the bearing and the shaft are situated at 
O and O’. The radial displacement OO’ is repre- 
sented by «. The shaft is moving to the left along 
a spiral orbit. The resultant of all the forces due 
to gas pressure is represented by the force vector 
Fx. The instant center of the orbit of O’ is situated 
at O’’. The line 0’ O”’ bisects the angle between 
Fx and the radius 0’ 0. 

It will be assumed that, over a restricted range, 
the force Fy varies directly as the radius 0 0’, and, 
hence, behaves like a spring. The force Fx may be 
resolved into three components. 


Fx cos- 6 Y(O 0’) 


which pulls the shaft toward the bearing center, 
Fx sin 6 


which accelerates the shaft center 


along the orbit, 
and a component 


Fx cos 6 sin 6 


which provides the energy required to stretch the 
“spring” 

The angular velocity w of O’ about 0 is given by 
the equation 


maw-(O O0’)=Y(0 O’)= Ye 


7 

\ Mn 

is the effective mass of the shaft in the 
plane of the bearing. 

The linear velocity w(0’’ O’) of the center of the 

shaft may be resolved into a tangential component 


where m 


a0" 0’ 


cos 6=— we 
and a radial component 


w(O0’’ O’) sin @=we tan #6. 

Thus the line of centers may be regarded as rotating 
about the center O with an angular velocity w, while 
the radial distance ¢ expands as it rotates. The rate 
of expansion is given by the equation 


where € is the value of the eccentricity when t=0. 
Thus if the direction of the resultant force lies ahead 


Figure 3. Whirl: decreasing amplitude. 


of the line of centers (@ positive), the diagram ex- 
pands exponentially with time; conversely if the 
direction of the resultant force lies behind the line 
of centers (6 negative) as shown in figure 3, the dia- 
gram shrinks in toward the origin and the bearing 
should be stable, even though the force Fx may be 
quite small. 

In practice, the forces due to the film of gas may 
be far from linear, and at a certain speed of the 
shaft, if the eccentricity is increased from zero, the 
sign of the angle @ may change from positive to nega- 
tive at some critical value of the eccentricity. Under 
these conditions the diagram should expand asymp- 
totically toward the critical value of the eccentricity 
and whirling of stable amplitude should occur. This 
is often observed in practice. 


3. Fluent Film Lubrication 


The action of a simple sleeve bearing is illustrated 


in figure 4. A shaft of radius r is surrounded by a 
sleeve of radius r+6, the interspace being filled with 
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Figure 4. Hydrodynamic pressure distribution. 





a continuous film of lubricant. The center of the 
shaft 0’ is displaced from the center O of the sleeve 
by a small distance e. If the shaft rotates counter- 
clockwise, the film of lubricant is carried into a con- 
verging channel on the upper right-hand side of the 
diagram, and out of a diverging channel on the lower 
left-hand side of the diagram. Thus the pressure is 
increased on the upper right-hand side and reduced 
on the lower left-hand side. The resultant force due 
to the pressure distribution is at right angles to the 
eccentricity. This is true, however, only for very 
small values of the eccentricity where compressibil- 
ity effects can be neglected; at larger values of the 
eccentricity, compressibility effects introduce an 
asymmetry and a component of force directed 
toward the center of the bearing is generated. 
Suppose the eccentric displacement of the shaft 
from the center of the bearing is maintained at a 
fixed value by means of a short link free to rotate 
about 0. The forces exerted by the film of gas on 
the shaft will cause the shaft to whirl about the 
center O. The whirling velocity w will be the veloc- 
ity at which the gas film is carried into the conver- 
> where w, is the angular velocity 


gent channel, i.e., 


of rotation of the shaft. 
Thus 


ie 
5" (1) 

This phenomenon which is known as “‘half-speed 
whirl” is often encountered in the bearings of high- 
speed spindles, particularly under conditions of light 
loading. 

If the shaft is of mass m, then there will be a 
centrifugal force mw*e due to the whirling which 
acts in the direction of the eccentricity. The 
pressure in the film of lubricant can only exert 
forces at right angles to this direction, and if half- 
speed whirl is occurring, this pressure is zero; 
consequently, if the short link is removed, the shaft 
will move outward along a spiral orbit until it comes 
into contact with the sleeve. As it moves out, 
however, compressibility and other effects cause 
the pressure distribution to become asymmetrical. 
Under these conditions, the new pressure distribution 
provides a force which can be resolved into a tangen- 
tial component and an inwardly directed radial 
component. In practical bearings, small deviations 
from perfection, such as variations in the film 
thickness, may introduce such a radial component 
even when the shaft is at the center of the bearing. 
At moderate speeds this radial component may be 
sufficient to balance the centrifugal force due to 
whirl. 

As the speed is increased, however, the radial 
component of the force due to the film increases 
as the first power of the speed, while the centrifugal 
force due to whirl increases as the square of the speed. 
Thus it is evident that at some critical speed the 
centrifugal force should become equal to the radial 
restoring force, and at this speed the shaft should 


become unstable. 
served in practice. 

In a perfectly constructed bearing, however, it 
is evident that a whirl of small amplitude should 
commence even at a very low shaft speed. This 
prediction has been verified experimentally by 
Boeker and Sternlicht [1956]. 

The tangential force Fy, which is generated in a 
bearing of radius 7, of radial clearance 6 and of 
length Z, is given by the equation 


This phenomenon is often ob- 


F 12rur’ Le 2 ( w 
“ & (2+n?)(1—n?)!2 
where 


is the viscosity of the gas in lb-sec/in.? 
; the eccentricity in inches 
; the eccentricity ratio 
; the rotational velocity 
; the angular velocity of whirl and 


L , 
f (<n is the end leakage factor which depends on 
#<5 


L 
— and n. 
r 


The first, second, and fourth terms on the right- 
hand side of eq (2) are taken from a paper by Ford, 
Harris, and Pantall [1957]. The third term, accord- 
ing to Barwell [1956], is introduced in order to 
account for the effects of whirl. Equation (2) 
holds only where compressibility effects can be ne- 
glected so that the pressure distribution is symmetrical 
about the line of centers. 

The value of the end leakage factor has been de- 
termined experimentally by Ford, Harris, and 
Pantall [1957]. According to the curves in figure 18 
of their paper, the value of the end leakage factor 
in bearings whose length is less than twice the di- 
ameter is approximately 


ABs tA . 
f( > ee Tr. \) 


and for small values of n, the leakage 
independent of the eccentricity ratio. 

Thus, at small values of the eccentricity ratio, the 
tangential force in a sleeve bearing whose length is 
less than twice the diameter is given approximately 
by the equation 


factor is 


Fy = Eye = 


Qrur* Le w ) 
> =i 
6° 


where Ey is the tangential “‘spring” modulus. 
Equation (4) may be written 


ae 

. ? rpr’ L?w,€ 2w 

Fy= Eye =~— ( 1— ) 
6 w, 


The force Fy can be either positive or negative 
depending on the value of w. In fact, Fy=0 when 
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™ 

is exactly equal to > 
on the shaft due to the gas film under this condition. 
The force Fy acts at right angles to the displacement 


There is no pressure acting 


WwW ’ ° ° ° . 
e, and when w<->> then Fy, acts in the direction of 


shaft rotation and whirl is accelerated; conversely, 


w ’ . : . . 
when w >>> then Fy, acts in the opposite direction 


to shaft rotation and whirl is retarded. 
An appreciation of the magnitude of the tangential 


force in a simple sleeve bearing at an eccentricity 
may be obtained from the following example. 


E’ram ple 1 


Radius of shaft 0.437 in. 

Radial clearance 50 X10"* in. 
Length of bearing 0.281 in. 
Viscosity of gas 2.97 «10 
Angular velocity of shaft 1.88 > 
Angular velocity of whirl Zero 


’ Ib-see/in.? 
10* radians/s 


3.14 2.97 10~"& (0.437)? x (0 


[Lox a os" 


Bil) 1268 > 


Eu 


Ea 784 lb in, (7) 


Thus at an eccentricity of 0.000393 in. (0.001 em) 
the tangential force is approximately 0.308 lb weight. 

It is evident that in order to insure stability at 
finite speeds of the shaft, the bearing must be de- 
signed so that a radial restoring (centering) force 
Fp is generated. if this force is proportional to the 
eccentricity, it behaves like a spring, and it may be 
expressed as 


Ye Fp (S) 


where Y is the modulus of the “spring” and ¢ is the 
eccentricity. An appreciation of the radial spring 
modulus required in one of the bearings of the 
miniature turbine may be obtained from the follow- 
ing example. 


Exram ple 2 


The effective mass of the shaft in the plane of the 
most heavily loaded bearing is 0.775 lb, and the 
designed rotational speed is 3,000 rps. The corre- 
sponding half-speed whirl valonke is 1,500 rps=9430 
radians/s. Inserting these values in the equation 
Ye mw’e we get 


Y =0.775 X 9430?/32.3 x 1; (9) 


178.200 Ib/in. (10) 


At an eccentricity of 0.00039 in., the centrifugal 
force would be approximately 69.5 lb weight. 


4. External Pressurization 
The necessary restoring force may be generated 
by means of external pressurization. The sleeve is 
provided with a number of equally spaced injector 
holes, as shown in figure 5. Compressed gas is 
supplied to the injector holes and flows radially out- 
ward in the clearance space between the shaft and 
the sleeve. The shaft is found to float centrally in 
the sleeve, a well-known effect. 

The reason for this effect becomes evident from an 
examination of a plane circular bearing as shown in 
the insert of figure 6. In this figure the relationship 
between the mean pressure between the plates and 
the plate separation is shown for various supply 
pressures. The graphs in figure 6 show the variation 
of pressure with the separation of the plates for 
plates of * in. diameter. In order to find the 
corresponding variation of pressure for plates of 
any diameter, the scaling relationship h?~D is used 
in conjunction with a constant orifice to pad diameter 
ratio. Accordingly, the graphs shown in figure 7 
are drawn for pads * in. diameter and a plate 
separation of 0.00165 in., which corresponds to a 
plate separation of 0.0015 in. with plates %. in. 
diameter. 

As shown in figure 5, the pads in this bearing are 
not circular but have a rectangular configuration, as 
shown. It will be assumed that each pad is equiva- 
lent to a circular pad %» in. diameter. The mean 
pressure in the clearance space between a pair of 
diametrically opposed pads [Sixsmith, Wilson, Bir- 
mingham, 1961] is shown in figure 7. Curves a and 
b show the variation of mean pressure with clearance 
on the right-hand and left-hand sides of the shaft. 
Curve ¢ shows the variation of differential pressure 
with the eccentric displacement of the shaft. Over 
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FIGURE Conventional gas-bearing arrangement. 
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Figure 6. Load versus plate separation curve. 
a restricted range the differential pressure (curve ¢) 
may be assumed to be linear, thus the corresponding 
restoring force is also linear, and is directly propor- 
tional to the displacement. 

If the direction of the radial displacement makes 
an angle y with a diameter through the center of a 
pair of diametrically opposed pads, the component 
of the deflection with respect to the pads is € cos y. 
Similarly, the component of the force generated by 
the pair of diametrically opposed pads in the direc- 
tion of the displacement is also proportional to cos y. 
Thus the total restoring force along any radius is 
given by 


1o0P 
2 Oe 


a) 7) 0 
A,e> *‘{cos? y|8, 


oP 
€ 
4 ~ O€ 


A, (11) 


where Y is the “spring”? modulus of the gas film, N 


is the number of pads, A, is the area of a pad, and 
1 oP 
2 O« 


as shown in figure 7. 


is the differential pressure due to a single pad 


> 


Exam ple v0. 


In the present bearing there are eight pads, and 
for this particular number of pads, 


>> [cos? y|]$,=4 


— 


(12) 
and (11) becomes 
oP 


Ye«=2A,;, >. (13) 
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FiGguRE 7. Mean pressure on gas-bearing pads. 


The effective area A, of a pad is 0.076 in.’, and the 
slope of curve ¢ at 200 psig supply pressure is 44,000 
lb/in.?) in. Putting these values in (13) we obtain 
the value of Y. Thus 


2 0.076 * 44,000 


6,700 Ib/in. (14) 
At an eccentric displacement of 0.00039 in., the 
radial restoring force is 2.61 lb. This is less than 
the corresponding centrifugal force (example 2) by 
a factor of 27. Thus it is evident that the radial 
restoring force, although quite large, is insufficient 
to prevent the shaft from spiralling outward until 
it touches the bearing. It is then necessary to devise 
a method whereby the centrifugal force can be 
reduced to a value below the available radial restor- 
ing force in order to obtain a stable bearing. This 
can be accomplished by reducing the shaft whirl 
speed. The method used to reduce whirl speed is 
presented in sec. 6. 


5. Analysis of Turbine Rotor Instabilty 


In an early design of a turbine [Sixsmith, 1959a] 
in which the turbine was mounted on a flexible shaft, 
it was discovered that a hydrodynamic exciting- 
force can be generated by the turbine rotor. At a 
certain critical pressure of the gas supply, a whirling 
of the shaft would commence. The shaft could be 
made to whirl even when it was prevented from 
rotating. In order to suppress the whirl it was 
necessary to provide damping in the form of oil-filled 
dash-pots. 
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Figure 8. Force distribution on turbine rotor. 


An explanation of this phenomenon is as follows. 
The effect appears to be due to forces which are 
generated by the variation of blade-tip leakage with 
rotor eccentricity. Consider the turbine rotor shown 
in figure 8. If the rotor is displaced downwards as 
shown, the tip leakage is reduced around the lower 
half of the rotor and correspondingly increased 
around the upper half. Around the lower half of the 
rotor the tip leakage is less, hence an increased 
thrust is developed on the blades. This force is to 
the right for counterclockwise rotation. Around the 
top half, tip leakage increases and the thrust on the 
blades is reduced. This can also be regarded as an 
added component of force to the right. Conse- 
quently there is a net force to the right acting on the 
rotor as a whole. This force is at right angles to 
the displacement, and therefore, it is in the same 
direction as the hydrodynamic force which is gen- 
erated in a simple sleeve bearing. 

An estimate of its magnitude may be made if we 
assume that the loss of torque due to tip leakage is 
directly proportional to the tip clearance, i.e., 


dT — [ _ 2arshe, ] 
27 WV(1{—T2) 


element of torque acting on 
between ¢ and ¢+d¢ 

total torque acting on the rotor 

tip radius—inches 

root radius—inches 

mean tip clearance between ¢ and 2+d¢ 
inches 

axial component of velocity of 
the tip clearance in./s 

axial component of velocity of 


(15) 


where 


the blades 


in. lb 


gas through 


gas through 
the annulus—in./s between r, and 7». 

If the rotor is deflected downwards a distance e 
the tip clearance h at a point which makes an angle 
with the direction of the displacement is given by 

h=6-4 


€ COS @¢. (16) 


The horizontal force dF to the right which acts on 
the small segment of blading lying between @ and 
o¢+d¢ is given by 


dF = cos @ 


Thus 


, Tdo 2arjhv 
dF:=5 1——.—— | cos ¢. 
221 m T(r —T13)V2 


° 


Substituting (6+ 


where 7,,=mean radius of rotor blades. 


e cos ¢) for h we have 


271m 


- (6+ € cos 6) | COS ¢. 


On integrating between O and 2z, all terms will 
vanish except the one containing cos’? ¢ thus we 
hive 

2ry?; i ie co 
21% _(} sin ¢ cos¢-+i 4) 
(71 —72)2 \2 Aiea 
Ter\v, l l 
; ; Qo 
TIm(11—1o)Vo\, 2 4 
Ter,v, 
‘m(Ti—T3) V2 


which may be written as 


where 


Thus it appears that the whirl-inducing force is 
proportional to the product of the torque and the 
eccentricity, but is independent of tip clearance as 
long as the assumption that the blade passage 
efficiency is inversely proportional to tip clearance is 
valid. Since for a given eccentricity it is propor- 
tional to the torque, it is capable of inducing a 
whirl even when the shaft is prevented from rotating. 

In the miniature turbine the value of the torque 
at 3000 rps and 3 kw of power is 1.406 lb-in. The 
values of 7, 7», and r. are 0.436, 0.397, and 0.36 in. 
Since the blade outlet angle is 30°, the value of 
es ; 
>, 8 2, approximately. 


Putting these values in (17) 
we have 
1.406 X 0.4386 K 2X € 
0.397 (0.436°—0.3607) 


F=51.15e lb. (18) 


Thus at an eccentricity of 0.00039 in. the tangential 
force is 0.02 lb weight. 
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the cavity, being in phase with the pressure, is like- 
wise lagging with respect to the current and eccen- 
tricity vectors. This force can be resolved into 
components whose directions are radial and tan- 
gential with respect to the eccentricity vector. 
The tangential component acts in the opposite 
direction to the hydrodynamic force due to the 
bearing and the turbine rotor. By suitable design 
of the bearing, the tangential component of force 
can be made to exceed the exciting forces so that 
the net tangential force is opposed to the tangential 
motion of whirl, and consequently, the whirling 
motion should be damped out. 

The radial component acts in the opposite direc- 
tion to the centering force due to the injector holes, 
and therefore it reduces the load carrying capacity 
of the bearing. This is inherent in the design and 
cannot be avoided. 

As shown in figure 9, a shaft of radius 7 is sur- 
rounded by a sleeve of internal radius r+ 6. Gas from 
eight injector holes each of radius 7;, equally dis- 
placed around the circumference of the sleeve, flows 
into eight cavities each of capacity C. The ma- 
jority of the gas escapes from each cavity through a 
leak orifice of resistance R. The remainder escapes 
at the ends of the bearing, but it will be assumed that 
this portion is negligible. Each cavity and leak 


Figure 9. Gas bearing with stabilizing cavities. 


6. The Analysis of a Method of Stabilizing 


a Gas-Lubricated Bearing 


In order to ensure stability, it is necessary to reduce 
the angular velocity of whirl until the centrifugal 
force becomes less than the radial restoring force. 
A reduction in the angular velocity of whirl may be 
effected if the bearing is designed to generate a force 
which acts in the opposite direction to the hydro- 
dynamic exciting forces due to the bearing and the 
turbine rotor. 

A bearing designed to produce this stabilizing 
force is shown in figure 9. In addition to the injector 
holes, as shown in figure 5, the shaft is surrounded 
by a number of pneumatic phase shift networks 
which are spaced midway between the injector holes. 
Each network consists of a cavity which acts as a 
pneumatic capacitor and a leak orifice which acts 
as a pneumatic resistor. A stream of gas flows into 
each cavity from the adjacent injector holes, and 
escapes from the cavity via the leak orifice. The 
pressure in the cavity is determined by the flow rate 
and the resistance of the orifice. In the event of a 
whirling of the shaft, the flow of gas into each cavity 
is modulated, the modulation being in phase with 
the eccentricity vector. As a consequence of the 
modulated inflow, the pressure in the cavity is like- 
wise modulated, but by virtue of the finite volume 
of the cavity, the phase of the pressure is lagging 
with respect to the phase of the current and the Figures 10-15 (inc.) Vector diagrams illustrating the action 
eccentricity. The force on the area of shaft facing of a stabilizing cavity. 
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behaves like an electrical resistance-capacity net- 
work: the capacity of a cavity being defined by the 
equation 
OP ,,0Op 
Ot ot 


p OP 


Pan 


(19) 


and the resistance of a leak by 


P 


R i, lb/s 


(20) 


where V isthe volume of the cavity, pis the density of 
the was in the cay ity, P is the pressure in the cavity, 
J. is the current of gas accumulating in the cavity, 
and J, is the current flowing through the leak: ¢ is 
time 

A cavity receives its gas supply from both of the 
adjacent injector holes. Suppose now that the 
shaft is whirling with an eccentricity ¢€ and a counter- 
clockwise angular velocity w. The motion of the 
center of the shaft may be represented by a complex 
vector ¢« rotatine counterclockwise with the 
angular velocity (fig. 10). 

Referring again to figure 9, it will be seen that 
the current of gas J entering cavity 1 will be given by 


same 


21) 
where vis the velocity of the gas at the edges of the 
injector holes and p, is its density. The values of 
the clearance at the injector holes adjacent to cavity 
l are given by 


€ COS ( wt 


-€ COS ( wt + 


where N is the number of cavities. 
The time ¢ is assumed to be zero when ¢ is pointing 
vertically downward through the center of cavity 1. 


Equation (21) may now be written 


£)+seooe (a+ )} 


q ‘ T 
‘p 4 26+-2e cos wt Cos ‘i . 
rn 


\ 


The current flowing into cavity 1 thus consists of 


a steady component J, which is given by 


I 27 


9,6 (23) 


and which rotates in 


vector « and whose 


an alternating component / 
phase with the eccentricity 


instantaneous value is given by 


IT, cos wt =2rr vp,e COS wt COS a (24) 

In the above analysis, the effects of shaft rotation 
and whirl have been neglected. These effects influ- 
ence the flow rate into the cavity. Referring to 
figure 9, it is seen that the clearance between the 
shaft and the bearing on the right-hand side of the 
cavity is greater than the clearance on the left-hand 
side. As a consequence of the rotation of the shaft, 
more gas is carried out of the cavity on the right- 
hand side than is carried in on the left-hand side. 
Thus there is a net flow of gas out of the cavity. An 
estimate of the value of this component of current 
may be made if it is assumed that the principle of 
superposition may be applied; i.e., that the velocity 
of the gas at any point in the film is the vector sum 
of the radial velocity from the injector holes and the 
tangential velocity due to hydrodynamic action. 

The instantaneous value of the hydrodynamic 
component of current J, is given by 


& sin wt L rpo(hs h 


w ‘§ e 

~ mdash | 5+ 
= 

| é +€ COS ( wt 


| t sin — 
w, Lr poe Sin wt sin => 
: A 


where J, is the alternating flow of current, w, is the 
angular speed of the shaft, is the radius of the shaft, 
L. is the effective length of the cavity, and p, is the 
density of the gas in the clearance space. 

Referring again to figure 9, it 
consequence 


is seen that as a 
of the whirling motion, the shaft is 
approaching the cavity and therefore the volume of 
the cavity is decreasing. This may be regarded as a 
current J, flowing into the ecavitv, the volume of 
the cavity being regarded as constant. The in- 
stantaneous value of the current is given by the 
equation 

ly sin wt 


2whrpoe sin wt sin y (26) 


where J, is the alternating flow of current and w is 
the angular velocity of whirl. 

Thus the resultant current J, due to the effects of 
shaft rotation and whirl is given by the equation 


° 2w : . Tv on 
[Sin wt Lr prew ( 1— ) sin wf sin V (27) 

This current may be represented by the vector Iq 
(fig. 11) which rotates 90° ahead of the eccentricity 


vector e. 
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The resultant current flowing into the cavity is 
represented by the vector I’. It is seen that 
I’=T, sec B. (28) 
The direct current J, flowing into the cavity may 
be represented by vertical line OA, figure 11, and 
the instantaneous value of the alternating current 
by the vertical line OB. This is the projection on a 
diameter through cavity 1 of the alternating current 
I’ which rotates ahead of the eccentricity vector 
with a phase angle £. 
The value of 8 is given by 


Te Lrpw 2w ra 
, ( 1— ) tan => 
I, 2xr,vp, w N 


tan 6 


(29) 


The steady current J, flowing through the resist- 
ance ? generates a pressure given by 


P,=]R. (30) 


This is represented by the vertical line OD, figure 12. 
The alternating current generates an alternating 
pressure component 
I R sec Bp 
\ 1+w*( ad , \ l 


LR sec B 
+-tan? (a+ £8) 


P’=['Z=—I 2 sec B 


(31) 


where Z is the impedance of the resistance-capacity 
network. Withrespect to the alternating component 
of current, the capacity and resistance are effectively 
in parallel and, consequently, the phase of the 
pressure lags behind that of the current and the 
eccentricity. The phase angle @ is given by 
tan (a+ 8B) =oRC. (32) 
The currents due to the alternating pressure P’ 
are represented vectorially in figure 13. The vector 
Ip represents the alternating current flowing in the 
resistance f, and the vector Ic represents the alter- 
nating current flowing in the capacity. The values 
of Ip and I¢ are given by the equations 


I / 


- 8 


P! eC=P' wl 


Pp 

P 
The resultant of these two currents must be equal 
to and in phase with the alternating current J’ 
flowing into the cavity. 

The pressures Py) and P’ acting on the effective’ 
area <A, of the shaft, which forms one wall of the 
cavity, generate forces 


(35) 


and 

P=: AF”. (36) 
These forces may be represented by the vertical 
vector OG (fig. 14) and the alternating vector OJ, 
respectively. The alternating vector F’ may be 
resolved into a component F’ cos @ which is in 
phase with the eccentricity and a component F’ sin 
a which is lagging 90° behind the eccentricity (fig. 15). 

These forces are generated by a single cavity and 
act along the radius from the center of the shaft 
through the center of that cavity. With a sym- 
metrical array of cavities, the resultant of the steady 
forces Fy will be zero; this, however, is not true of the 
alternating components. 

The resultant force on the shaft at any instant is 
found by adding vectorially the instantaneous values 
of the forces due to each of the cavities. 

The instantaneous value of the force due to cavity 
l is 


F’ cos (wt—a). (37) 


This is the projection of the alternating vector F’ 
along a diameter through the center of the cavity. 
At a given instant it is represented in magnitude and 
direction by the vertical line ON as shown in figure 
i4. 

The polar diagram of the force due to each of the 
eight cavities is shown in figure 16. At an instant 
of time t, the vectors numbered 1 to 8 represent the 
forces due to the corresponding cavities. 

The polar diagram represents a field of force ro- 
tating with the angular velocity of whirl. This is 
illustrated in figure 17, which shows the values of 
the eight forces due to the cavities when ¢ has in- 
creased to t+ At. These forces can be combined into 
a resultant which by svmmetry will lie along OJ. 
The magnitude of this force may be written as Fs 





FIG.17 


FiGuRES 16 Polar diagrams 


stabilizing cavities. 


of fore generate d by 


Experimental testing has shown A/2 to be 
value of A» as illustrated in figure 20. 


i conservative estimate for, the 
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and is given by 


F.=- (38) 


Substituting from 


becomes 


(36), (31), (29), and (24) this 


N rAgrjvp,Ze cos v 


sec 6 


(39) 


Fp x@Ze (40) 


where 


ui 


Ne sec B. 


G=NrAz jUP; COS 


(41) 


G is a constant which depends on the design of the 
bearing and the speed of the shaft. 
vector GZe may be resolved into components @GZe 
cos a in phase with the eccentricity vector or GZe 
lagging 90° behind the eccentricity vector. 
The component of force GZe cos a due to the stabi- 
lizing cavities is in the same direction as the centrif- 
ugal force due to whirl and in the opposite direction 
to the centering force due to the injector holes 
defined in (11 Consequently, the centering force 
is reduced to Ye—GZe cos a. Whirl can set in at 
an angular velocity given by 


Sill @ 


Mmawe Ye GZe COS Q@; 


(42) 


Y —GZ cos a\'? 
w—( ) 


(43) 


Mi 


These forces are represented as vectors in figure 18. 


ENTER OF BEARING 


RE 1S. Illustration of forces on shaft. 


The alternating 


The component GZe sin a provides a force at 
right angles to the eccentricity vector and acts in 
the opposite direction to the forces due to hydro- 
dynamic action or excitation in the bearing and the 
turbine rotor F (17). This is the damping force 
which is required to provide stability. Provided it 
is greater than the sum of the bearing and turbine 
rotor excitations, the direction of the resultant of 
all the forces lies behind the line of centers, and any 
whirl, such as might be excited by an external 
shock, is attenuated until its amplitude becomes 
zero. 

An estimate of the magnitude of the hydrodynamic 
force due to the bearing may be made as follows. 
The bearing is regarded as a simple sleeve bearing 
in which the shaft rotates with an angular velocity 
w, and a whirl velocity w. The value of the whirl 
velocity w is given by (43). The hydrodynamic 
force which is generated at an eccentricity € can be 
calculated with the aid of (5). The force thus 
calculated may be expressed in the form 

Eye=F y (44) 
for small values of the eccentricity ratio «/6. Ey is 
a constant whose value depends on the design of the 
bearing. In the above calculation no account is 
taken of the disturbance in the flow pattern due to 
the various holes and slots in the surface of the 
bearing. These have the effect of reducing the 
pressures due to hydrodynamic action and therefore 
the calculated value of the hydrodynamic force Eye 
is an overestimate. Thus it errs on the safe side. 
The force Eye is represented vectorially in figure Is. 

The rotor excitation kT is given by (17). This 
applies however to the force in plane of the rotor; 
hence, the value must be multiplied by some factor 
k,, depending on the overhang, in order to give the 
force in the plane of the bearing. This force may 
thus be expressed in the form 

F,=k' Te (45) 
where 

k’ =k, k. (46) 
This force is shown vectorially in figure 18. The 
term k’T can be regarded as the rotor excitation 
reduced to the plane of the bearing. Referring to 
figure 19, it is seen that the excitation due to the 
rotor is greater in the inner bearing than it is in 
the outer bearing. For the inner bearing the value 


of k; is given by 
L; 
k,=( > } (47) 
L. T Lg 
The stabilizing force GZe sin a which is generated 
by the bearing must be greater than the sum of the 
two excitations outlined above if the shaft is to be 
stable. 
The condition required for dynamic stability is 
thus given by 
GZe sin a>k’ Te 


Eve ° 


(48) 
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Figure 19. Shaft resonant system. 

Provided this condition is satisfied, the magnitude of 
e will progressively decrease to zero and the vector 
diagram in figure 18 will shrink into the origin. 
On the other hand, if the sum of the excitation is 
greater than the damping term, the vector diagram 
will progressively expand until « becomes equal to 
the radial clearance—at which point the bearing 
may be destroyed. 

In the above analysis the conditions which deter- 
mine the dynamic stability of the bearing have been 
outlined. However, the bearing must also possess 
static stability, i.e., it must be capable of supporting 
a radial load, such as the weight of the shaft if it 
were horizontally supported. In this case, w becomes 
zero and if w, is also zero then @ is zero, and since tan 
(a+ 8)=wCR, the phase angle a also takes the value 
of zero. The stabilizing force @Ze sin a thus reduces 
to zero, and the decentering force GZe cos a reduces 
to the value GRe. The centering force now becomes 

F,=(Y—GR)e, (49) 
and it is evident that GR should be less than Y in 
order to ensure static stability. If GR should be 
greater than Y, the shaft will be pushed to one side 
of the bearing. It is evident that the value of the 
leak resistance must be carefully chosen; it must be 
large enough to ensure dynamic stability; if it is 
made too large, however, the shaft will not float. 

If we write 

’p 

B, (50) 
then the static stability condition will be satisfied if 
B is less than one. The static centering force may 
now be written 


F,=Ye(1—B). (51) 


The value of the leak resistance is given by 


. Dee 
R —* 
I Ge 
Substituting for G from (41) we get 


BY 


NrAgr vp; cos 


< sec B 


In the present bearing, B was chosen to be 0.75, 
and consequently, the static load carrying capacity 
is 0.25 Ye. The area of the shaft A, which faces each 
cavity is estimated to be as shown in figure 20. The 
value of 8 is given by (29). Asa first approximation 
it may be assumed that B=0 and sec B=1. The 
flow velocity v is assumed to be sonic at the throat 
formed by the annulus 2776; ie., the discharge 
coefficient is assumed to be unity. The diameter of 
the leak orifices may be calculated with the aid of 
this equation. The discharge coefficient for a 
circular orifice is approximately 0.6. 

Since the shaft is supported on two bearings, each 
of which behaves like a spring, it has two principal 
modes and frequencies of whirl. The calculation of 
the modes and frequencies of such a resonant system 
is described in standard textbooks on mechanical 
vibrations. However, for the convenience of the 
reader, the procedure in this specific application is 
outlined below. 

The shaft of the miniature turbine and its supports 
are shown schematically in figure 19, where the 
following symbols and the corresponding values are 
used. 


Mass of shaft—lb mass 

Radius of gyration of shaft about the center of 
gravity—in. 

Stiffness of gas film in outer bearing 

Stiffness of gas film in inner bearing 

Distance of outer bearing from 
eravity—in. 

Distance of 
eravity—in. 


Ib/in. 
Ib/in. 
center of 
inner from 


bearing center of 


It is assumed that the two radial spring constants 
are equal, i.e., 


ie Y; "— GZ cos a. (54) 


In this equation, a value for a which should 
provide the optimum damping is required. This 
value is given by the equation 


90°—B 
c 9 


a= 


where 


Tn 




















Figure 20. Effective area facing cavity. 





as shown in figure 11. But J, cannot be evaluated 
until the value of w is known, and this depends on 
the value of cos a. 

This difficulty may be circumvented if it is as- 
sumed as a first approximation that the angle 6 is 
zero and the angle a@ is 45 With these values of 
a and 8, Z becomes equal to R 2 and (54) becomes 


(57) 


Y— GR. 


In the present calculation the 


arbitrarily chosen so that 


value of PR is 


(58) 


The two frequencies of oscillation are given by 
the equations 


Latotalt @—ort(®) 


where 


mr; 
Equations (59) and (60) are approximate and do 
not include the gyroscopic effects associated with 
rotation of the shaft. However, such effects can be 
ignored for the case considered here. 
the nodes from 
gravity are given by the equations 


The distances of the center of 


. (62) 
a 
positions of the two nodes are illustrated in 
19. From the positions of the nodes it is 
at the predominant amplitudes of whirl in 
r and outer bearings occur at the lower and 
frequencies of whirl. 

values of @, and w, may now be employed 
ind (54) in order to find more precise values 
nd },. From the revised values of Y, and 
more refined calculation of the values of a, 
nav be made. However, in view of the many 


uncertainties in the computation, the first approxi- 
mation is usually sufficiently accurate for practical 
purposes. 

The value of the phase angle 8 may now be 
calculated with the aid of (29). The optimum value 
of the phase angle a is given by (55). 

The values of the stabilizing cavities may now be 
calculated. We have 


wVpR 


tan (a+ 8) P 


wR (63) 
(a B)P 
w ph 


V, tan (64) 
+-8)P 
wopR 


V, tan (a 


(65) 


where w, and w, are the whirl angular velocities of 
the outer and inner bearings, V; and V, are the 
corresponding volumes of the stabilizing cavities, and 
p is the density of the gas at P units of pressure. 

Referring to figures 12 and 11, it is seen that the 
vectors of pressure and current are given by the 

equations 
Y’=['E (66) 

and 
ie | 


sec B. (67) 


The value of Z is given by the equation 


R 
ttan” (a+ 8) 


Z (68) 


The value of the force F’, figure 
single cavity is given by the equation 


14, due to a 


Pp’ =A?’ (69) 
where A, is the effective area of the shaft 
which forms part of the 
value of the force GZe 
equation 


(fig. 
walls of the cavity. 
(fg. 18) is given by 


20) 
The 
the 
GZe— F’. (70) 


Thus the value of @GZe is given by 


NARI sec Bp 


GZ 
te 9 1+ tan? (a+) 


(71) 
The corresponding tangential and radial 
ponents of force are GZe sin a and GZe cos a. 


The condition for dynamic stability is given by 
(48); Le., 


com- 


GZe sin a >k’ Te+ Eye. (72) 
Provided this condition is satisfied, the direction of 
the resultant of all the forces acting on the shaft 
lies behind the line of centers as shown in figure 18. 
Thus the orbit of the center of the shaft should be 
a spiral of decreasing amplitude as shown in figure 3. 
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FIGURE 21. diagram of 


rrocedure, 
l 


) , 
forces illustrating design 


Consequently, the bearing should be stable. This 
was proved experimentally in the bearings of the 
miniature turbine which is described by Birmingham, 
Sixsmith, Wilson |1962] 

The design procedure and the relationship between 
the various forces may be illustrated with the aid 
of the geometrical construction shown in figure 21. 
The angle 8 is preferably chosen to be less than 
20°, assuming zero whirl velocity, and the diameter 
of the injector holes is calculated according to (29) 
to satisfy this condition. 

In figure 21 we have: 


OS=BYe 


O’D=BYe sec B 


AN=1/2BYe sec B 


O'H=AN—O’'A sin B 


1/2BYe see BCI (76) 


sin B). 


Thus the stabilizing force O’H may be expressed 


directly in terms of the centering force Ye. The 
condition for dynamic stability may now be expressed 
as 
1/2BYe see B(1—sin B) >k’ Te4+- Eve. (77) 
In a trial design the number of pads and the 
radial clearance are adjusted until the above sta- 
bilitv condition is satisfied. Then the diameter of 
the leak orifices and the volume of the stabilizing 
cavities may be calculated as follows. 


720—591—64 4 


The values of Y, and Y, according to (57) are 
given more simply by the equation 


(78) 


and then the values of the whirl angular velocities 
may be calculated with the aid of (59) and (60). 
The diameter of the leak orifices may be calculated 
according to (53). Then the values of the stabilizing 
cavities may be calculated according to (64) and 
(65). This completes the design. 

In the theory outlined above, it is assumed that 
the leakage along the shaft from the stabilizing 
cavities is negligible. Such leakage introduces a 
phase shift which tends to promote instability. 
Conversely, a leakage along the shaft and into the 
stabilizing cavities tends to promote stability. 
Consequently, it is desirable that the pressure at 
the ends of the bearing should be equal to or greater 
than the static pressure in the stabilizing cavities. 


7. Conclusion 


It is possible and quite practical to provide stable, 
externally pressurized gas-lubricated bearings for 
high-speed shaft operation. The stabilizing cavities 
employed in this design are a definite asset to a 
bearing system of this type and have been shown 
by actual tests to suppress shaft whirl inherent in 
externally pressurized bearings. 

The bearing described here is adaptable to small 
shafts (one inch diameter and smaller) which are 
lightly loaded and require speeds well beyond the 
allowable operating range of rolling contact bearings. 
For miniature applications, the size of the bearing 
and shaft system is limited only by practical fabri- 
cation techniques. 

This type of bearing is particularly advantageous 
for miniature turbine expanders used in cryogenic 
refrigerator systems, where the process gas can be 
used as the bearing gas medium. In this arrange- 
ment, the problem of contaminating the process gas 
from the bearing lubricant is eliminated. It also 
eliminates a mechanical seal between the bearings 
and the process expansion system. 

A prototype turbine expander, incorporating ex- 
ternally pressurized gas-lubricated bearings, was 
fabricated and tested with helium gas as the operat- 
ing medium. The bearings were designed using the 
criterion described in this manuscript. The tests 
showed that the bearings were reliable and performed 
quite satisfactorily. No detectable shaft whirl was 
encountered up to the design speed of the turbine 
shaft. 

8. Summary of Nomenclature 
A,=Area of pad—in.’? 
A,=Effective area of shaft 
cavity—in.? 
Constant defined in eq (50) 


surface facing 





Capacity of cavity 

Diameter of pad—in. 

Tangential ‘spring’? modulus—lb/in. 

Horizontal force on blading—lb 

Static centering foree—lb 

Force due to gas pressure on shaft—lb 

Tangential foree—Ilb 

Radial restoring foree—lb 

Steady force—lb 

Alternating vector of forcee—lb 

Vector sum of forces due to N cavities—lb 

Constant defined in eq (41) 

Mean tip clearance between ¢ and ¢+d¢ 
in. 

and bearing. 


between shaft 


See figure 9—in. 


( ‘learances 


Alternating component of current flowing 
in cavity—lb/s 

Current through leak—lb/s 

Flow of gas into cavity—lb/s 

Steady component of flow into cavity 
lb/s 

Alternating component of flow into cavy- 
ity lb/s 

Alternating flow due to shaft rotation (eq 
ne 
aw 

Alternating flow due to shaft whirl 

Alternating flow due to the combined effects 
of shaft rotation and whirl 

Alternating flow into cavity 

Constant defined in eq (17) 

Correction for overhang of rotor 

Constant defined in eq (46 

Length of bearing—in. 

Effective length of stabilizing cavity—in. 


see figure 19 


Mass of shaft in plane of bearing—lb 
Number of pads in bearing 


€ 


Eecentricity ratio 

6 

Center of bearing 

Center of shaft (see fig. 1) 

Instant center of orbit of shaft (see fig. 2) 

Pressure in cavity—psi 

Steady component of pressure in cavity- 
psi 

Alternating 

cavity 


component of pressure in 


Resistance of leak orifice—sec/in.? 

Radius of shaft—in. 

Turbine rotor tip radius—in. 

Turbine rotor root radius—in. 

Turbine rotor mean radius—in. 

Radius of gyration of shaft 
center of gravity—in. 

Total torque acting on turbine rotor 

time 


about its 


in.-lb 


sec 


- Volume of stabilizing cavity (see eq (19))— 
in.® 
Volume of stabilizing 
bearing—in. 
Volume of stabilizing 
bearing—in.° 
Velocity of gas at edge of injector hole—fps 
Axial component of velocity of gas through 
rotor tip clearance—fps 
Axial component of velocity of gas through 
the annulus between 7; and r,—fps 
Radial spring modulus—lb/in. 
“Spring” modulus of film in 
bearing—Ilb/in. 
“Spring” modulus of 
bearing—lb/in. 
Z=Impedance of resistance capacitance 
work 


cavities in outer 


cavities in inner 


outer 


Fas 


as film in inner 


net- 


Phase angle between complex pressure 
vector and eccentricity € 

Phase angle between resultant 
and eccentricity 

Angle between ¢ and a diameter through the 
center of a pair of diametrically opposed 
pads 

Radial clearance—in. 

Shaft eccentricity—in. 

See figure 2 

Viscosity of gas—lb-sec/in? 

Density of gas in cavity—lb/ft 

Density of gas at edge of injector hole— 
lb/ft 

Density of gas in clearance space between 
shaft and bearinge—lb/ft 

See figure 8 

Angvlar velocity of whirl—radians/s 

Angular velocity of whirl in outer bearing— 
radians/s 

Angular velocity of whirl in inner bearing— 
radians/s 

Angular speed of shaft 


eurrent J’ 


radians/s 
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Infrared Spectra of Asphalts: Some Aspects of the 


Changes Caused by Photooxidation 
Paul G. Campbell and James R. Wright 


(November 29, 1963) 


The infrared spectra of coating-grade asphalts to which model organic compounds 


were added were 
asphalts. 


compared with spectra of the original asphalts and with photooxidized 
By this method, it was shown that the principal oxygen-containing products 
formed in the asphelt during air- or photooxidation were aldehydes, 
The presence of esters could not be confirmed. 


ketones, and acids. 


The effects of air-blowing conditions and elevated temperature storage on several as- 


phalts were examined by infrared spectroscopy. 


of air-oxidized apshalts, 


Elevated temperature (~260 C) storage 


under nitrogen, decreased absorbance in the carbonyl band at 


1700 em~! and reduced durability to accelerated weathering. 


Photooxidation reactions that occurred at the surface 


of thick (625 w) asphalt films, as 


detected by infrared spectroscopy, were the same as those which took place in the thin 


25 uw) films. 


Also, the photooxidation rates were comparable. 


Only negligible amounts 


of photooxidation occurred at a depth of 10 uw beneath the surface. 


1. Introduction 


The initial application of infrared spectroscopy 
to the analysis of aphaltene-type materials was 
made a decade ago by Hadzi [1] ' who studied quali- 
tatively the infrared spectra of the light-petroleum- 
solvent-insoluble, benzene-soluble fraction of various 
materials. Friedel and Queiser [2] made an infrared 
spectrum of an apshalt film in 1956 and used it to 
illustrate the similarities to the spectra of coal 
asphaltenes. Also, in 1956, Knotnerus [3] applied 
infrared analysis to determine the oxygen-con- 
taining functional groups formed during the air 
blowing of asphaltic bitumens. 

Stewart [4] described the infrared spectra of as- 
phalt fractions separated by chromatographic 
techniques in 1957. In this research, asphalts 
were separated into four fractions and an infrared 
spectrum was made of each fraction. Samples of 
the original asphalts were then exposed for 30 
months outdoors, subsequently fractionated, and 
an infrared spectrum was made from each fraction 
of the exposed asphalt. Some of the structural 
changes produced by weathering were elucidated 
by this method. Stewart worked with solvents, a 
heated cell, potassium iodide pellets, and mineral 
oil mulls in obtaining the spectra of various asphalts 
and asphalt fractions. 

In a study of the absorbance of asphaltic materials, 
Schweyer [5], in 1958, made infrared spectra of asphalt 
films 0.1 mm thick. For wavelengths below 8 uy, it 
was necessary to use carbon tetrachloride solutions 
of the asphalts. Schweyer attempted to explain cer- 
tain differences in the properties of asphaltic mate- 
rials on the basis of the ultraviolet and infrared 


Figures in brackets indicate the literature references at the end of this paper 


spectra. Romberg et al. [6], used infrared spectros- 
copy in the analysis of asphalt components from gel, 
sol-gel, and sol asphalts in 1959. From the results 
thus obtained, the percentages of various carbon 
linkages in each fraction were calculated and com- 
pared with other properties such as C/H ratio, 
elemental analyses, and molecular weight. 

Beitchman [7], in 1959, described a method for the 
preparation, irradiation, and infrared analysis of 
thin films of air-blown asphalts of about 25 u thick- 
Beitchman observed that increases in absorb- 
ance at 343! 1 (2.91 nw), 1700 em™ (5.88 uw), and 
1030 em (9.71 yw), (OH, C=O, C—O groups, re- 
spectively], were produced by exposure of unsup- 
ported films of asphalt to the radiant energy of a 
carbon arc. Yen [8], in 1962, published the results 
of an extensive investigation into the structure of 
petroleum asphaltenes by infrared analysis. 

The research described by Beitchman [7] was the 
first successful effort to produce infrared spectra of 
asphalt films without resorting to solvents, mulls, 
pellets, or salt-crystal supports. The unsupported 
film was particularly adaptable to research on asphalt 
photodegradation. The films could be scanned, irra- 
diated, and rescanned without modification to the 
asphalt other than that caused by photodegradation. 
Wright, Campbell, and coworkers modified the un- 
supported film technique [9-11] and used it exten- 
sively to determine photooxidation rates of asphalts 
[10], the effect of temperature and relative humidity 
on asphalt oxidation [11], and the relationship be- 
tween carbon-are intensity and asphalt oxidation 
[12]. Changes in absorbance at 1700 em, the car- 
bonyl index, were employed to measure photooxida- 
tion. 

The present paper is an extension of this research 
with the general objective of interpreting more fully 
some of the changes observed in the infrared spectra 
of asphalts due to photooxidation. Specific areas 


ness. 
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investigated include the nature of the carbonyl com- 
with absorbance at 1700 em spectral 
produced by air oxidation and by elevated 
temperature storage of asphalts, oxidation rates of 
thick films, and the depth of the photooxidation 
reaction in asphalt. 


peunds 


oOhanges 


2. Infrared Spectra of Asphalts Containing 
Model Organic Compounds 


Infrared 

yesence 1m a Mm iterial 
hydroxyl 

cos 

on the 

<ture, the 


cives evidence of the 
of characteristic groups suc h 
‘arbonyl groups. If the material i 
such information is often defini- 
other hand, the material is a 
nterpretation Is more complicated 
the location of any functional groups, relative to any 
pound in the mixture, is not normally known. 
Che complexity of : isphalts limits the interpreta- 
f their infrared spectra. Similarly, the changes 
ich occur In the infrared spectrum of an asphalt 
photooxidation are subject to limited interpreta- 
For example, with 
carbonyl band at 17¢ 
e spectra of 
considerable 
rbonvl group, 
id anhvdride, 


1 °2] 
Led]. 


spect I" ISCOP\ 


Con 


to absorbance in 
relative changes 
can be measured 
the source of the 
aldehyde, ketone, 
not be established 


respect 
0 cm, 
+} 


asphalt films 


witt ACCULAC) 10]: 
whether it be 
ester, etc., Can 
possible, however, 
identity of 
which 


npounds to an 


to gain some insight as to 
compounds by an indirect 
adding model organic 
scanning a film of the 
the changes produced in the 
ared spectrum. New peaks or 
uulders are produced if the compound type differs 
the asphalt Existing absorbance 
peaks are reinforced if the model-compound additive 
corresponds to compounds of the same types already 
present. in the asphalt. The model-compound- 
additive method was emploved in this research to 
the types of compounds formed during 
“idation of asphalt with special emphasis on 
nges which oceut 1700 


these 
consists of 
asphalt, 
ure, and obser, ing 


pproach 


absorbance 


those in 


ucidate 


Oto 


in the absorbance at 
2.1. Experimental Procedure 
a. Methods 


irradiation, and 
asphalt films has 
The model organic 
incorporated into the asphalt in 
ilowing manner: a known amount of a selected 
lt was heated on a hot plate under a blanket 
until melting occurred. The organic 
und was added to 2% by weight) and the 
it mixture stirred until relatively homo- 
The mixture was then covered and cooled 
temperature. The whole operation required 


‘he method for the preparation, 
red analysis of unsupported 
described elsewhere [10]. 

were 


pounds 


ns (25 uw thick) of the additive-treated asphalts 


were prepared aus described in {10}. Infrared scans 


were made of these films. For direct comparison, 
the spectrum of the additive-containing asphalt was 
made on the chart paper which contained the 
spectrum of the corresponding untreated asphalt 
from a film of similar thickness. Perkin-Elmer 
infrared spectrophotometers (Models 221 and 137) 
were used in making the infrared scans. 


b. Materials 


Two air-blown roofing asphalts were used in the 
present study: One asphalt (No. 16) was of Mid- 
continent United States origin and had shown a 
durability to accelerated weathering (ASTM D529- 
59T) of 93 days using the 51-9C cycle (9 min of 
cold water spray in each hour). The other asphalt 
(No. 9) was a California Coastal asphalt fluxed with 
reclaimed lubricating oil; its durability was 32 days. 
Properties of these asphalts have been described 
elsewhere [9]. 


2.2. Infrared Spectra of Air-Blown Asphalts 


a. Aspha!ts With No Compounds Added 


The infrared spectra of a typical air-blown asphalt 
before and after photooxid: ition are shown in figure 
1. Significant changes in absorbance may be 
observed at 3436 em 1700 em~', and 1600 em. 
According to Stewart [4], respective absorbances at 
these wave numbers are due to: OH. stretching, 
C=O, aromatic and conjugated alkene and OH 
bonding. Also, absorbance changes at other wave 
numbers have been reported by Beitchman [7] 
and Stewart [4]. 

Of the mentioned above, relative 
changes due to photooxidation are best measured at 
1700 em™ since this, the carbonyl band, is relatively 
free from interference by absorbances of other 
groups [14] and carbonyl absorbances can be deter- 
mined with greater accuracy than other groups such 
as OH and C—O. This band was used extensively 
by the authors in previous investigations on the 
photodegradation of asphalts [9-12], but the identifi- 
cation of the types of organic molecules responsible 
for the absorbance at 1700 em 


absorbances 


was not made. 


WAVENUMBER 


4000 3000 2000 


ABSORBANCE 


MIDCONTINENT 
ASPHALT 


— BEFORE EXPOSURE 


TER EXPOSURE 


WAVELENGTH , 4 


FIGURE | Infrared spectra of a midcontinent U.S.A 


before 


. asphalt 


and after carbon-arc trradiation. 
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b. Asphalts With Organic Compounds Added 


The increased carbonyl absorbance observed in 
asphalt upon photooxidation may be due to any 
carbonyl-containing organic compounds. To_ in- 
vestigate this phenomenon by the method described 
in section 2a, the asphalts were mixed with esters, 
aldehydes, ketones, acids, and anhydrides, respec- 
tively. A total of 14 different compounds was used. 
These compounds are listed in table 1 together with 
their chemical structure, characteristic absorbances, 
and the absorbances imparted to asphalts containing 
“ach compound. 


TABLE 1. 


Charac- 
teristic 
ibsorb- 


ance (s) 


Compoun 


cm cm 
1730-1717 | 1730-1720 
1310-1250 1275 (broad) 
1150-1100 | 1115 (broad 


1750-1735 
1172 


"1717 


1300-1100 broad 


broad) 


shoulder 


2700-2500 
*1708 
1350-1180 


*1704 
1300-1260 


Benzaldeh ye 


1 Source L 
the compound. Unmarked entries refer to the class of compounds. 


Model organic compounds added to asphalts No. 


Absorbance (s) 
in asphalt 2 


shoulder 


broad) 


broad) 


J. Bellamy, The infrared spectra of compler molecules, 2d ed. (John Wiley and Sons, New York, N.Y., 1958 


The effects of various organic compounds on the 
infrared spectra of asphalt No. 16? (with and without 
organic compounds added) are shown in figures 2 to 
7. The changes produced by phthalic anhydride 
are shown in figure 2. Two carbonyl absorbance 
bands, characteristic of anhydrides are present at 
about 1845 and 1775 em™!. Bellamy [15] reports 
values of 1845 and 1775 cm“! for phthalic anhydride. 
In addition to the carbonyl absorbance, there is evi- 
dence of the C—O—C stretching vibration at 1250 
em. According to Colthup [16], vibration in the 
range of 1310-1210 cm™ is characteristic of cyclic 
anhydrides. 


9 and No. 16 and their contribution to the infrared spectra 


Structure 


Compound 


cm 
1705-1685 
1440-1325 


975-82 


1705-1685 
1230-1160 


n-Heptalde- 2 1740-1720 
hyde 1440-1325 


Acetophenone *1686 


*1690 
1325-1215 


*1848 
Maleie anhy- *1790 


lrice ‘ 1310-1210 broad) 


*1845 
Phthalic anhy- *1775 
j 1310-1210 


dride 


Asterisk (*) indicates value is specific for 


* Reference is made to the regions of increased absorbance in the asphalt-additive mixture over that found in the original asphalt. 


Asphalt 16 was selected for this study because it has a high initial absorbance 
in the carbonyl! band and would reflect only large contributions due to additives. 
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WAVENUMBER , CM Similarly, citral, an  a-8-unsaturated aldehyde, 
increased carbonyl absorbance at 1705 em™ (fig. 4). 
a-8-Unsaturated aldehydes normally absorb in the 
range of 1705-1685 em™! [18]. 

The addition of a fatty acid, stearic acid, to asphalt 


3 x 200 hel 
AL - = a 

\ SM” MARA 16 resulted in spectral changes typical of this class 
\] \ 7 ‘(hw Ge eee of compounds [19]. Slightly increased absorbance 
t| a oa | was observed in the 2700-2500 em™ range (bonded 
if — ASPHALT NO. 16 OH stretching vibration), at 1705 em~! (C O 
- BLEND: ASPHALT 16- vibration), and in the broad region around 1250 
The increased absorbance in the 1705 
>is as «Oem! region enhances that of the original asphalt. 
The addition of esters to asphalt 16 produced an 
effect not previously observed with any of the other 
classes of carbonyl compounds in air-blown asphalts 
Changes in the infrared spectrum of asphalt No.16 or in photooxidized asphalts. Two absorbance 
used by the addition of phthalic anhydride. peaks were formed due to C=O vibration. This 
effect is shown graphically in figures 6 and 7 with 
n-butyl sebacate and tristearin, a_ triglyceride, 

respectively. 

Characteristics C=O absorbance for alkyl esters is 
in the range of 1750-1735 em™ with the n-butyrates 
at the lower end of the range [17]. This absorbance 
range is notably higher than the carbonyl frequency 
for alkyl ketones and aldehydes due to the influence 
of the adjoining oxygen [14]. The C=O absorbance 
for n-butyl sebacate in asphalt 16 is 1740 em™ 
while the absorbance of the asphalt alone in this 
region is at 1700 cm™!. Similarly, the respective 
absorbances for C=O in tristearin and asphalt 16 
are 1750 em ~ and 1700 em~!. The triglycerides 
absorb in the range 1751-1748 em™! [21]. 

The C—O stretching vibrations due to esters 

. : (1300-1100 em~') are found in the asphalt spectra 
nee op 1 atapol teea a an No. 16 with both n-butyl sebacate and tristearin. With 
ee ee the sebacate, an absorbance is found at 1185-1175 

em7!. Sebacates absorb at 1172 em™'! [20]. Charac- 
‘iii is teristic C—O stretching modes for triglycerides are 
1250, 1163, and 1110 em [21]. With tristearin, broad 
1000 90¢ : absorbance occurred at 1250, 1160, and 1110 em™!. 
— Photooxidized asphalts also exhibit increased ab- 
sorbance in the 1300-1100 em™! region (see fig. 1). 
However, with photooxidized asphalts, none of the 
characteristic bands for C—O vibrations in esters can 
be distinguished. 





$+ 





; 
+4 
} 


| 
1} : PHTHALIC ANHYDRIDE em”! (fig. 5). 


—— ASPHALT 0.16 . ° 
SE a 2.3. Discussion 
CITRAL 

5 — By adding model organic compounds to an asphalt 
; and observing the changes produced in the infrared 
spectra, it was shown that ketones, aldehydes, and 
acids reinforced the absorbance assigned to the C=O 
Changes in the infrared spectrum of asphalt No. 16 frequency at 1700 cm ', while acid anhydrides and 
caused by the addition of citral. esters did not contribute in this wave number region. 
The anhydrides exhibited two carbonyl absorbance 
bands at 1845 and 1775 cm™!. The esters developed 
iddition of mesityl oxide, an a-6-unsaturated an absorbance in the 1750-1735 cm~! region forming 
to asphalt 16 produced increased absorbance a doublet with the carbonyl absorbance in the 
at 0 em and directly reinforced the normal asphalt. Thus, it appears that the increased absorb- 
absorbance of the asphalt at this wave number ance found in oxidized asphalts is due largely to 
(fig Thompson and Torkington [17] show the ketones, aldehydes, and acids with no contribution 
frequency of the stretching vibration of the carbonyl from anhydrides and little or no contribution from 
group mesityl oxide to be about 1690 em™!. esters. This generalization is applicable to asphalts 
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oxidized in air, in oxygen, and in oxidants other than 
oxygen [22] as well as specimens subjected to photo- 
oxidation. 

The insignificant contribution of esters as a source 
of oxygen in oxidized asphalts is worthy of further 
comment. In 1955, Goppel and Knotnerus [23] 
published a paper on the fundamentals of bitumen 
blowing in which more than 60 percent of the oxygen 
in the blown bitumen was attributed to ester groups.* 
Furthermore, since esters link two different molecules, 
it was postulated that esters contributed to the 
formation of higher molecular weight materials. 
Proof of ester formation, based on saponification 
values and spectroscopic evidence, was described in 
a separate paper by Knotnerus [3]. These two papers 
have been referenced frequently in the subsequent 
literature dealing with asphalt oxidation under a 
variety of experimental conditions [24-27]. 

Since the present work does not support the 
contention that esters are the predominant oxygen- 
containing molecules in oxidized asphalt, the dis- 
crepancy may be explained as follows: The spectro- 
scopic evidence offered by Knotnerus [3] is based 
upon the nearly complete disappearance of a strong 
carbonyl absorption band at 1724 cm™! and the 
formation of a new band at 1587 em™! (carboxylic 
acid salts) when blown bitumens are saponified with 
excess alkali. (The remaining absorbance at 1724 
cm~! was believed due to aldehydes and ketones.) 
The possibility of high molecular weight acids was 
ruled out because the hydrogen bridges were not 
observed between the solvent employed in the spec- 
troscopic method and the potential acids. This 
evidence, coupled with saponification values, was 
offered as proof that esters were the predominant 
molecular type. It is the authors’ contention that 
the presence of ‘‘acid-solvent hydrogen bridges” 
would be extremely difficult to detect by infrared 
spectroscopy in a solution containing the whole 
asphalt. And furthermore, that the absence of such 
bridges under these conditions is rather weak evi- 
dence for selecting esters rather than acids as the 
predominant oxygen-containing molecule. 

Saponification values also may be misleading be- 
cause compounds other than esters may produce 
positive result. Siggia [28], in discussing the limita- 
tions of the saponification equivalent method, pointed 
out that the presence of reactive compounds would 
give erroneous results to the method used in ester 
determination. According to Shriner et al., [29], 
hot concentrated alkali affects functional groups 
other than esters. For example, aldehydes with 
a-hydrogen atoms undergo aldol condensation, alde- 
hydes with no a-hydrogen atoms undergo Canniz- 
zaro’s reaction, and 8-diketones undergo cleavage 
under the influence of hot alkali. 

To supply further evidence of this fact, saponifica- 
tion equivalents were measured on n-heptaldehyde 
and acetophonone for 4 and 24 hr alkali treatment. 
The values for the aldehyde were 525 and 253, and 
for the ketone were 10,020 and 1,469, respectively. 


3 Similar effects were postulated for aged bitumens, 
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Figure 5. Changes in the infrared spectrum of asphalt No. 16 


caused by the addition of stearic acid. 
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RE 6. Changes in the infrared spectrum of asphalt No. 16 


caused by the addition of n-butyl sebacate. 
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The present work does not, of course, rule out ester 
formation entirely. It is well known that a-6-un- 
saturated esters exhibit carbonyl absorbance in the 
region of 1720 em™!. The addition of ethyl cinna- 
mate (C=O frequency of 1717 cm™') to asphalts 
16 and 9 broadened the carbonyl absorbance but 
gave no doublet as was the case with normal satu- 
rated esters. Also, in photooxidation of asphalts, 
poorly defined shoulders are sometimes observed on 
the main carbonyl absorbance region; these might be 
attributed to ester formation since asphalts contain 
both acids and aleohols. It should be noted that 
esters containing one or more 6-hydrogen atoms in the 
alkyl group may decompose photochemically by an 
intramolecular rearrangement into an olefin and the 
corresponding acid [30]. In general, it is concluded 
that the oxygen in oxidized asphalts is principally in 
the form of aldehydes, ketones, or acids (the present 
results do not distinguish among these classes) with a 
relatively minor part as esters. 


3. Spectral Changes Caused By Air-Oxida- 
tion and Elevated-Temperature Storage 


The changes in the infrared spectra of the asphalts 
described in section 2 were produced by the addition 
of model organic compounds to the asphalts. Most 
hard asphalts, however, are manufactured by air 
oxidation at temperatures in the order of 260 °C 
(500 °F). Therefore, it was of interest to elucidate 
the spectral changes produced in asphalt under these 
conditions as well as those produced by storage of 
the asphalt at the elevated temperatures used in the 
air-oxidation process. 

Three asphalt fluxes, of Midcontinent United 
States, Venezuelan, and Southeast United States 
origin, were used. The designations of these fluxes, 
physical properties, the conditions of air 
oxidation, and the softening points (S.P.) (ASTM 
D36-26) of the oxidized asphalts are given in table 2. 
(For more details on these materials, see Greenfeld 
[31] 


some 


Infrared spectrum of Kansas L 


ce asphalt flux, 
softening point, 39° 


Film supported on N 


The infrared spectrum of the Kansas L-22 asphalt 
flux, before air oxidation, is shown in figure 8.1. The 
area around 1700 cm7! shows practically no absorb- 
ance indicating little or no oxygen in the C=O form. 
Upon air oxidation to a softening point of 104 °C 
(219 °F), considerable absorbance was observed in 
the 1700 em™! region indicating the formation of 
C—O oxygen (see fig. 9). 

Absorbance around 1700 cm~! may be attributed 
to C=O in the form of aldehyde; ketone, and/or 
acid. The change in absorbance at 1300 em7! may 
well represent alkyl ketones and aromatic aldehydes. 

Another change appeared in the region ° of 760— 
720 em. Well-defined absorbances occurred at 
750 cm~! and 725 em™ in the asphalt flux. In the 
air-oxidized asphalt, these absorbances were less 
intense and poorly defined. 

The infrared spectrum (not shown) of the air- 
oxidized but softer asphalt (S.P.=80 °C) was similar 
to that of the more highly oxidized asphalt (S.P.-= 
104 °C). The main difference was one of absorb- 
ance intensity in the 1700 cm™! region; the softer 
product showed less absorbance and, consequently, 
less carbonyl formation. 

Storage of the In-grade ° 


Kansas L-22 asphalt 
(S.P. 


104 °C) at 243 °C (470 °F) for 19 hrs under 
a nitrogen atmosphere caused a decrease in 5.P. to 
95 °C. Correspondingly, the absorbance intensity 
in the 1700 em~! region decreased, but no shift in the 
absorbance region was observed. Apparently some 
thermal degradation of the oxygen-containing mole- 
cules took place. The storage temperature of 
243 °C was severe enough to cause decarboxylation 
of organic acids. 


ABSORBANCE 


Ficure 9. Infrared spectrum of coating-grade asphalt (soften- 


ing point, 104° C. 


Made from Kansas L-22 flux 


‘ The softening point of the L-22 flux was 38 °C (101 °F 
ported film. For this reason, the 
chloride plates 

Bellamy [32] interprets strong absorbances in this region as due to aromatic 
s" bstitution Yen and Erdman [8] assign absorbances for asphaltenes in the 
865-760 cm~! range as due to aromatic out-of-plane bending frequencies. 

The designation *‘In-grade”’ refers to an asphalt which meets the industry’s 
specifications with respect to softening point and penetration. 


; too low for an unsup- 
isphalt film was prepared between sodium 
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The relationship between asphalts modified by 
air-oxidation or elevated-temperature storage and 
their photooxidation rates was determined. The 
physical properties of the asphalt fluxes before and 
after asphalt oxidation are shown in table 2. The 
Physical 


TABLE 2. properties 


after 


of asphalt fluxes before and 


air oxidation 


Before air 
oxidation 


After air oxidation Air-oxidation 


conditions 
Designation 

Soft- 

ening 

point 


Pene- 
tration 
at 25°C 


Below 
grace 


In- After Tem- 


Air 
grade storage perature Time 


rate 


1/10 


mm Softening point 


( 
38 249 
4) 152 


40) 190 


Kansas L-22 
Tia Juana 612 


46 
Talco 16-RF 


243 


194 


Conditions given are those used to produce asphalt designated “ 
Stored for 19 hr at 243 °C, under nitrogen 

Stored for 21 hr at 246 °C, under nitrogen 
‘Stored for 21.5 hr at 243 °C, under nitrogen 


In grade’’ 


experimental method was that employed by Wright 
and Campbell [10] in which 25 u films of each asphalt 
were exposed to the radiant energy of the carbon 
are. Oxidation was measured as a function of 
exposure by the changes occurring in C=O ab- 
sorbance at 1700 em~'. The asphalts exposed and 
the resultant carbonyl indices are shown in table 3. 

Oxidation-rate curves are plotted for Kansas 
L.-22 asphalts in figure 10. The asphalt that was 
least susceptible to photooxidation was the In- 
grade asphalt (S.P==104 °C). Film failure oc- 
curred after 16-hr exposure. The Below-grade 
asphalt (S.P.—80 °C) exhibited a slightly higher 
photooxidation rate, but film failure did not occur 
until after 18.5 hr. The least durable asphalt, 
with a much higher photooxidation rate and a 
film life of only 12 hr, was the asphalt subjected 


TABLE 3. 


Carbony] index 


Asphalt Kansas L-22 


Softening point, ~¢ 


104 
E xposu 


0.014 
027 
031 


O72 


135 
165 


Irradiated Atl 
Irradiated for 12.5 hr 
Irradiated for 18.5 hr 
‘One sample only 
Irradiated for 7.75 hr 
Irradiated for 17.25 hr. 


ir 
in 


is SMC-R Weather Ometer at 49°C and 40 percent R.H. for time 


0.18 


KANSAS L-22 ASPHALT 

0 - IN GRADE (S.P. 104 °C) 

4 - BELOW GRADE (S.P 80 °C) 

O- STORED 19 HR AT 243 °C 
(S.P 95 °C) 


O.16+ 


i2 
TIME , HOURS 


FIGURE 
RH. 
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10. Effect of carbon-are exposure 19 °C., 40 percent 


on the change in infrared absorbance of various Kansas 
22 asphalts at 1700 em-'. 


to elevated temperature storage. It appears that 
degradation reactions initiated thermally might 
well have continued under conditions for photo- 
degradation. 


4. Spectral Changes in the Surface of Thick 
Asphalt Films Upon Photooxidation 


4.1. Spectral Changes 


The spectral changes in the asphalt described in 
sections 2 and 3 were determined with thin films of 
25 uw (1 mil) thickness. However, the usual methods 


for the exposure of asphalts to solar or carbon-are 


Changes produced in carbonyl indices of various asphalts by irradiation ' 
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114 
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040 
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0.014 
034 
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076 
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» periods indicated 





radiant energy employ asphalt films of 625 uw (25 mils) 
thickness [33]. In order to compare by infrared 
methods the changes observed in the thin films with 
those which occurred on and near the surface of the 
thick films upon photooxidation, slices of 25 u 
thickness were removed from the surface of irradiated 
thick films and their infrared spectra determined. 

A Southeast U.S.A. asphalt, No. 6, was used in 
this study. An asphalt specimen was prepared by 
pouring the molten asphalt into a small aluminum 
weighing dish to a depth of about 625 uw. A 3.8-em 
diam microtome specimen holder was then embedded 
into the soft asphalt. Upon cooling, the asphalt 
was cut to the diameter of the specimen holder and 
the surface was shaved smooth with a microtome 
(Spencer, Model 860). 

The asphalt samples were exposed to the radiant 
energy of a carbon are at 49 °C and 40 percent R.H 
for the desired time period. The samples were 
mounted in the same respective positions in the ac- 
celerated weathering machine as were previous sam- 
ples of thin films. After exposure, a 25 u slice was 
removed from the surface, weighed to check thick- 
ness, and then pressed between sheets of unlacquered 
cellophane in a heated press to prepare a specimen 
suitable for infrared analysis. 

The infrared spectrum of a 25 u film of asphalt 6, 
prepared from the irradiated surface of a 625 uw film 
is shown with the spectrum of an unexposed thin 
film asphalt 6 in figure 11. The changes in the 
infrared spectrum of the asphalt obtained from the 
surface of the thick film are the same as those pro- 
duced in the thin films, namely those at 3435 em7! 
(OH), 1700 em=! (C=Q), and 1030 em! (C—O). 
The carbonyl absorbance at 1700 em7! and the 
absence of absorbance at 1740 cm! indicates that 
the oxygen is in the form of acids, aldehydes, and 
ketones. The slight increase in absorbance at 2750 
cm! suggests the presence of carboxylic acids and 
that around 1300 em the presence of alkyl ketones 
and aromatic aldehydes. 

Figure 12 shows the changes which occurred in 
the infrared spectrum of a thin film of asphalt 6 
under the same exposure conditions as those used 
with the thick films. Significant changes in absorb- 
ance occurred at 3435 em7!, and 1700 em, and 
increases in general absorbance at 1030 em™!, as was 
the case with the thick films. Even the absorbance 
ntensities in these regions were about the same. 
These results indicate that the photooxidation re- 
actions occurring in the thin films, as determined by 
infrared spectroscopy, are essentially the same as 
which take place at the surface of thicker 
asphalt films under the exposure conditions used 
the present work. 


those 


4.2. Photooxidation Rates 


\ comparison was made of the photooxidation 
rates of thin films of asphalt 6 and that which 
occuired on the surface of thick films of the same 


asphalt. To do this, a series of 625 uw films were 
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Figure 11. Comparison of the infrared spectra of unexposed 
asphalt No. 6 (25 p-thin film) with that of a film made from 
the surface (25u-microtome slice) of the exposed asphalt (625p 
thick film). 

Exposure was for 7.5 


hr to carbon-are radiation at 49 °C and 40 percent R.H) 
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FicureE 12. Infrared spectra of asphalt No. 6 (25 wu thin film) 
before and after 7.5 hr. eX posure 


49°C and 40 percent R.H. 


to carbon-are radiation at 


prepared on microtome holders and exposed for 


varying durations to carbon-are radiant energy at 
49 °C and 40 percent R.H. After selected exposure 
periods, 25 slices were removed from each film and 
specimens prepared for infrared analysis as described 
in 4.1.' 

A survey of the exposure data (table 4) shows that 
no significant difference exists in the photooxidation 
rates of the two films based on carbonyl index for at 
least the first 10-hr exposure. For greater durations, 
the surface became hardened, brittle, and uneven. 
As a result, films of uniform thickness could not 
be sliced from the asphalt specimen. This was 
manifested by the large variance of AA values 
obtained beyond the 10-hr exposure period. For 
this reason, photooxidation rates of 25 w and 625 
films could not be compared for an exposure period 
of more than 12 hr. 


7In order to calculate the carbonyl] index (AA at 1700 cm-') [10], the initial car- 
bony! absorbance of these films was estimated from a graph (previously prepared 
on thin films of asphalt 6) of film thickness versus initia] carbonyl absorbance. 
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TABLE 4. Changes produced in carbonyl indices of 25 w films 
and in 625 w films of asphalt No. 6 by irradiation! 


Carbony! index 
Exposure 
time 
surface of 
625 yw film? 


25 yw film 


0.011 0.015 
024 022 
045 048 
064 059 
101 096 
116 4.125 


Irradiated in Atlas SMC-R Weather-Ometer at 49 °C and 40 percent R.H 
for time periods indicated 


2 Based on a 25 w microtome slice from surface of 625 w asphalt film. 
Irradiated for 7.5 hr 


4 Average of three value 


s, 0.148, 0.085, and 0.132. 


4.3. Depth of Photooxidation in Asphalt 


In a separate experiment employing the micro- 
tome, asphalt specimens were exposed to carbon-arce 
radiation for 10 hr at 49 °C and 40 percent R.H. 
A 10 uw slice was removed from the surface and then 
a second slice of 25 yw thickness removed. The 
initial carbonyl absorbance (1700 em!) of this 
second film, i.e., 10 to 35 uw below the exposed surface, 
was not significantly different from films of unex- 
posed asphalt 6 prepared in the usual manner [10]. 
The implication presented is that photooxidation 
of asphalt is apparently negligible at depths of 
>10 uw from the surface. This is in agreement with 
the conclusions of Dickinson et al., [34]. 


The authors express their appreciation to S. H. 


Greenfeld, Research Associate, Asphalt Roofing 
Industry Bureau, who kindly furnished the asphalt 
samples and physical property data on the asphalt 
fluxes. 
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Dielectric behavior of the flm formed on mica cleaved in 
moist air, S. Ruthberg and L. Frenkel, J. Res. NBS 68A 
(Phys. and\Chem.), No. 2, 173-183 (Mar.-Apr. 1984), 70 cents. 
Water absorbed on a freshly peeled mica cryst#l causes the 
loss tangent, D, to increase by 1 to 2 orders of magnitude. 
The nature of the film is investigated as a function of relative 
humidity by the measurement of D for the frequency range 
100 to 50,000 c/s with a capacitor geometry comprising con- 
centric, parallel electrodes of different diameter on opposite 
sides of the dielectric sheet This geometry is then analyzed 
as consisting of the two regions of that within the plates where 
E is normal to the dielectric plane and that at the edge wher« 
tangentisl FE exists. The first is considered in terms of an 
equis ilent eireult for a two lavered dielectric. The second 
is considered in terms of transmission line concepts. It is 
predicted and verified that the adsorbed film causes the first 
components to vary es 1/f and the second as yf where tis the 
thickness of the erystal. Numerical solutions are used to 
derive the behavior of A, R, and C of the adsorbed film itself. 
D for the i! direction is ~ 0.4 and follows a frequeney 
dependency like fresh snow Resistivities normal to and 
parallel to cleavage are considerably different while C is 
much less than expected for a surface film. It is suggested 
that the surface film is not continuous but instead is localized 
in patches, 


norm 


On the measurement of dielectric losses and surface con- 
ductivity of dielectrics in parallel plane test capacitors, L. 
Frenkel, J. Re VBS 68A Phys. and Chem.), No. 2, 185 

18S Va \ pr 196 70 cents. 

Thin slabs of dielectric materials are often tested for their 
dielectric properties In pl me pare llel plate condensers. When 
surface conductivity is present, as for instance in freshly split 
mica, losses not connected with the bulk of the material arise. 
The present paper deals with the general theory of such meas- 
urements. The system is reduced to an assembly of lumped 
elements superimposed on distributed transmission lines. 
The treatment includes the presence of possible airgaps under- 
neath the plates of the test condenser. It is shown that such 
losses depend on the reciprocal square root of the frequency. 
Losses due to this effect cannot be eliminated by guard ring 
messurement and much of the published data on the losses 
in mica must be reexamined in the light of the present work. 
Similar considerations may apply in the other 
materials. 


ease of 


Second spectrum of 
VBS 68A (Phys. 
1964), 70 cents. 

A preliminary 
published 


tungsten (Wi), D. D. Laun, J. Res. 
Chem.), No. 2, 207-252 (Mar. Apr. 


] 
and 


report on the second spectrum of tungsten, 
in 1958, presented 27 even energy levels and 50 
odd energy levels that were derived from 500 W_=11 lines 
ranging in wavelength from 1961.43 A to 4348.13A. The 
present paper submits data on 62 even levels, 132 odd levels, 
and 2,150 classified lines of W u, ranging in wavelength from 
1756.6 A to 6219.77 A. The ground state of the W* ion is 
represented by the ®D,4 level of a sextet D term arising from 
the 5d’‘6s! electron configuration, but the level intervals and 
magnetic splitting factors indicate considerable departure 
from LS-coupling, suggesting that 7j-coupling or intermediate 
coupling may be more appropriate for the W 11 spectrum. 


Optimal matchings, and degree-constrained subgraphs, A. J. 
Goldman, J. Res. NBS 68B (Math. and Math. Phys.), No. 1, 
27-29 (Jan.—Mar. 
The characterization of maximum-cardinality matchings in 
linear graphs, by the nonexistence of augmenting paths, has 
been extended by several authors to a similar characterization 
of maximum degree-constrained subgraphs. This paper 
contains a proof of the extended version by direct reduction 


1964), 75 cents. 


to the case of matchings. Possible algorithmic implications 
of the reduction are suggested. 


integral 
Math. 


Gaussian wave functions for polyatomic molecules: 
formulas, M. Krauss, J. Res. NBS 68B (Math. 
Phys.), Vo. 1, 35-41 (Jan.—Mar. 1964), 75 cents. 
Kxplicit formulas are given for the molecular coulomb in- 
tegrals that arise for all functions These 
results should expedite computational efforts with these basis 
functions. 
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and 


Gaussian basis 


Some basic microwave phase shift equations, KR. W. Beatty, 
Radio Sci. J. Res. NBS USNC-URSI 68D, No. 4, 345-354 
lpr. 1964), $1.00. 
The phase differences between terminal variables (voltage, 
current, or traveling wave types) at the output with respect 
to similar ones at the input of a 2-port are expressed in terms 
of the scattering coefficients of the 2-port and of the reflection 
coefficients of the system into which it is inserted. 
\ variety of phase shifts may be defined for a given 2-port, de- 
pending upon which of the terminal variables are considered, 
whether or not generator and load reflection coefficients are 
anish. and in which direction the 2-port has been 
inserted into the svstem. 
It is shown that a reasonable choice for one of the two “‘charac- 
teristic” shifts of a 2-port is Yo, the argument of S 
one of the seattering coefficients. It follows that the other 
“characteristic’’ phase shift is Jv The corresponding 
change in characteristic phase shift for variable phas« 
is the change in either Yy» or Yo (whichever is appropriate) 
from the initial to the final setting. 
Ideal phase shifters are discussed, and expressions for the 
change in output level of variable phase shifters are given. 
The importance of using a nonreflecting system in phase 
shift measurements is emphasized. 
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A simplified theory of diffraction at an interface separating 
two dielectrics, J. Kane and 8. N. Karp, Radio Sci. J. Res. 
VBS/USNC-U RSI 68D, No. 3, 303-310 (Mar. 1964), $1.00. 
Many electromagnetic problems involving more than one 
dielectric medium are not susceptible of an exact solution, 
when the appropriate boundary conditions are considered. 
The purpose of the present paper is to formulate a new 
boundary condition, which is capable of leading to mathemat- 
ically tractable problems, with limited sacrifices in accuracy. 


On the long term phase stability of the 19.8 ke/s signal trans- 
mitted from Hawaii, and received at Boulder, Colo., A. H. 
Brady, Radio Sci. J. Res. NBS/USNC-U RST 68D, No. 3, 
283-289 (Mar. 1964), $1.00. 

The use of VLF signals for intercontinental frequency com- 
parison has become very popular in recent years, and it has 
been shown by other workers that a precision of about 2 parts 
in 10 can be achieved with measurements over a 24 hr 
period. Phase records made at Boulder, Colo. of the NPM, 
Hawaii, 18.6 ke/s transmission have been studied for several 
periods of nearly two weeks duration in 1962. Deviations 
from an assumed linear frequency difference between the 
transmitter and receiver oscillators show an attainable pre- 
cision of 2.5 parts in 10!' in a 24 hr observation, extending to 
3.1 parts in 10” in a 192 hr observation. Without further 
data on the remaining differences between the oscillators, 
there is no evidence that the propagation conditions over the 
path are limiting this precision. It is pointed out that the 
50° (peak to peak) semiannual variation in the day to night 
change in phase, if attributed entirely to one level of reflection, 
would have an effect on precision of frequency comparison of 
about 1 part in 10”, 


Precise phase and amplitude measurements on VLF signals 
propagated through the Artic zone, F. H. Reder, C. J. Abom, 
and G. M. R. Winkler, Radio Sci. J. Res. NBS/USNC- 
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URSI 68D, No. 3, 275-281 (Mar. 1964), $1.00. 

The scope and purpose of an Arctic VLF propagation study 
are outlined, and results are presented for the period from 
February to June 1963. Typical diurnal phase and ampli- 
tude patterns of NPM and NPG transmissions received at 
Stockholm are given for month. The diurnal NPG 
amplitude change at Stockholm reached a maximum of 38 db 
during March. Monthly averages of NPM-Stockholm and 
NPG-Stoekholm diurnal phase changes are compared with 
those for other transmission paths. Twenty-four hour aver- 
ages of frequency differences of VLF transmissions were 
measured at various reception sites and intercompared. 
Residual errors ranged from 10 to 3X 10-", with standard 
deviations between 8.5 and 3.3 10 Arctic paths were 
found to be only a factor of about 2. Preliminary 
analysis of geomagnetic disturbance data and observed VLF 
phase anomalies indicates reasonable correlation for the 
NPG-Stockholm path. Finally, phase measurements on 
NBA taken simultaneously at Kiruna and Stockholm are 
compare d An SCN A event recorded by a Kiruna riometer 
clearly coincided with an NBA phase anomaly observed at 
Kiruna while the Stockholm NBA phase record showed no 
anomaly. 
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A light-modulated scattering technique for diffraction field 
measurements, A. M. Vural and ID. K. Cheng, Radio Sci. 
J. Res. NBS/USNC-U RSI 68D, No. 4, 355-362 (Apr. 1964 
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Propagation of electromagnetic waves through a continuously 
varying stratified anisotropic medium, G. H. Price, Radio Sci. 
J. Res. NBS/USNC-U RSI 68D, No. 4, 407-418 (Apr. 1964 
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Methods for the dynamic calibration of pressure transducers, 
J. LL. Sechweppe, L. C. Eichberger, ID. F. Muster, FE 
Michaels, and G. I VBS Mono. 67 (Dec. 12, 1963 


Paskusz 


his Monograph is designed to t the practicing engineer 
d with the problem of making 
brations of the rapidly ch 
vehicle technology 
der to calibrate and use dynamic pressure transducers, 
practicing engineer needs to be familiar with (1) the char- 
istic differential equations and their solutions, (2) the 
ids of analyzing pairs of input and output functions to 
mine the transfer function and the frequency response 
3) the methods of generating precise input functions, 
the specific methods for using precise experimental 
rements to determine the dynamic characteristics of 
ilar pressure transducer. 
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anging pressures in missiles and 


Compilation of the melting points of the metal oxides, S. J. 
Schneider, NBS Mono. 68 (Oct 10, 1963), 25 cents. 

\ compilation has 
metal 


been made of the melting points of 70 
xides published prior to January 1963. Both the 
nelting point and the equivalent value based on the 


ional Practical Temperature Scale of 1948 are pre- 


sented. Included in the survey is information on pertinent 
experimental details such as the method of temperature meas- 
urement, purity, furnace type, and environmental conditions. 


Weights and measures standards of the United States, a 
brief history, L. V. Judson, NBS Misc. Publ. 247 (Oct. 1963), 
(Supe rsedes Sci. Pap. 17 and Misc. Publ. 64.) 

A historical account is given of the standards of weights and 
measures of the United States from the time of the American 
Revolution through the year 1962. Current and historical 
standards of length and mass now in the possession of the 
National Bureau of Standards are listed and described. 


35 cents. 


Calibration and test services of the National Bureau of Stand- 
ards, NBS Misc. Publ. 250 (Nov. 22, 1963), 70 cents. 

This publication is a listing of the numerous calibration and 
testing services provided to science and industry by the Na- 
tional Bureau of Standards. The Bureau promotes accuracy 
and uniformity of measurement through its program of meas- 
urement services, including the calibration and testing of 
standards and standard instruments. An up-to-date listing 
of the Bureau’s calibration and test fee schedules was printed 
in recent issues of the Federal Register, with a large propor- 
tion of the fees being changed at that time. This publication 
contains all of this material, but in larger, more legible type, 
and fully indexed. Besides listing all NBS calibration serv- 
ices and their cost, the publication includes a statement of 
the Bureau's statutory functions, testing policy, and routine 
for securing the Bureau's calibration and test services 
original request to the reporting and use of test results. 


from 


Bibliography on ignition and spark-ignition systems, G. F. 
Blackburn, NBS Misc. Publ. 251 (Nov. 22, 1963) 15 cents. 
S ipersedes VB 

Approximately 30 references to books, papers, and reports 
are listed on ignition of combustible mixtures and 
ignition apparatus. The ignition of gases includes ignition 
by electric sparks and ares and by hot surfaces. The ref- 
erences on ignition apparatus are for the most part on ignition 
systems and components on internal-combustion 
with spark plugs listed separately from other 


S Circ. 580 
‘ 
I 


gaseous 


engines, 


components. 


Survey of the literature on safety of residential chimneys and 
fireplaces, H. Shoub, NBS Misc. Publ. 252 (Dec. 17, 1963) 
145 cents. 

\ search was made of recent technical literature on chimneys 
and fireplaces to determine what information was available 
on questions concerning their safety from fire hazard. The 
survey covered types of construction and materials proposed 
or currently in use with several types of fuel. The reports 
included investigations of heat transmission through chimney 
and fireplace walls, the effects of thermal shock on flue 
materials, and the required spacing of combustibles from the 
heated walls The bibliography presented brief 
reviews of some of the works 


? 


contains 


Recommended unit prefixes; Defined values and conversion 
factors; General physical constants, VBS Misc. Publ. 253 
Oct. 1963 wallet-si 
Wallet-size pocket card listing selected general physical con- 
stants recommended by the National Academy of Science- 
National Research Council and adopted by the National 
Bureau of Standards. Units are given in both the centimeter- 
gram-second System and the International System (MKSA). 
Reverse side lists unit prefixes recommended by International 
Committee on Weights and Measures and adopted by the 
National Bureau of Standards. Also lists defined values and 
conversion factors for several commonly used measurement 
units, 


5 cents e plastic card). 


Photometry of projectors at the National Bureau of Standards, 
L. Chernoff. NBS Tech. Note 198 Dec. 17, 1963), 25 cents. 
A projector consists of a light source and optical system which 
produces a beam of light. In the photometry of projectors, 
the intensity of the projector as a function of angle of viewing 
is measured. For many years the National Bureau of Stand- 
ards has conducted photometric tests on various types of 
projectors, and a variety of techniques and equipment has 
correspondingly been developed. 
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Photometric testing is carried out on a photometric range of 
which there are two at NBS: a short range on which the test 
distance may be varied up to a maximum of 30 meters and a 
longer range with a fixed test distance of 279 meters. Photo- 
metric measurements are made by comparison with a stand- 
ard lamp of known luminous intensity in a specified direction. 
These comparisons are made with photosensors which are 
color-corrected by filters so that the spectral response is 
similar to the CIE luminous efficiency function. While 
precision requirements are not severe, considerable care is 
required to keep the experimental errors within the desired 
limits. 


Carbon resistors for cryogenic liquid level measurement, 
R. C. Muhlenhaupt and P. Smelser, NBS Tech. Note 200 
(Oct. 9, 1963), 25 cents. 

Data are shown in graphical form. One set of plots presents 
resistance ratio Rg¢/R, as a function of “‘warming up” time at 
various levels of constant power dissipation. A second set of 
plots presents resistance ratio R,/o as a function of nominal 
resistance at various levels of constant power dissipation. 
The use of the data and the design of a practical liquid level 
indicator are discussed in the appendix. 


A technique for extrapolating the 1 ke values of secondary 
capacitance standards to higher frequencies, R. N. Jones, 
NBS Tech. Vote 201 (Nov. 5, 1963), 15 cents. 

A simple technique is described for extrapolating the 1 ke 
values of certain two-terminal capacitors to higher frequencies 
without incurring serious losses in accuracy. The method is 
intended for use with air capacitors having binding post or 
banana plug type connectors. Because of inherent errors, 
such connectors are not appropriate for measurement where 
the highest accuracy is required. For this reason, it is recom- 
mended that calibration of this type be performed by second- 
ary laboratories and not be submitted to the National Bureau 
of Standards. 


BOUMAC, A macro-programming system 
putation, J. H. Devenney and J. J. 
203 ( Dec. 18, 1963), 30 cents. 

The BOU MAC computing system will allow a scientist having 
no programming knowledge to perform the more commonly 


for scientific com- 
Sopka, NBS Tech. Note 


used macro computations such as matrix operations, solution 


of linear numerical 
statistical analyses, etc., using a high-speed electronic com- 
puter. The scientist needs only to know the various mathe- 
matical operations which are required to obtain the informa- 
tion he wants from the data and to learn to use a simple code 
with which to specify these mathematical operations. 


equations, integrations, curve fitting, 


The normal phase variations of _ 18 ke/s signals from NBA 
observed at Frankfurt, Germany, A. H. Brady, A. C. Murphy, 
and D. D. Crombie, NBS Tech. Note 206-1 (Dec. 8, 1963), 25 
cents. 

Observations of the normal phase variations of the 
signals radiated from the Canal Zone 
furt, Germany are given in the form of monthly averages and 
standard deviations at 5 minute intervals. The relation 
between mean diurnal phase variations and the diurnal varia- 
tion in path illumination is shown. The mean diurnal height 
change is 21 km which is reduced to 20 km during summer and 
winter and increased to 22 km during the equinoxes. Short 
term normal phase differences are also briefly discussed. 


IS ke/s 
and received in Frank- 


Solar heating, radiative cooling, and thermal movement— 
their effects on built-up roofing, W. C. Cullen, NBS Tech. 
Note 231 (Dec. 16, 1963), 25 cents. 

Twenty different built-up roof construction 
covered with five surfacing materials, were subjected to 
natural solar heating and nighttime cooling. The tempera- 
tures and temperature changes observed during winter and 
summer exposures are discussed. The data indicate that the 
temperature attained in a roof membrane is influenced by 
the absorptance and emissibity of the surface as well as the 
thermal and physical properties of the substrate to which 
the roofing is applied. The data show that roofings placed 
over insulation may be heated to a temperature of 80° F. 


specimens, 


above ambient due to solar heating and may be subcooled 
as much as 20° F. below the ambient due to radiative cooling. 
The thermal movements which occur in the components of 
a roof system due to temperature change are discussed in 
relation to built-up roof performance and failures. Thermal 
expansion data are presented for some composite bituminous 
membranes. The data show that these membranes undergo 
greater thermal movements than most other components of a 
roof system and the rate of expansion is not linear but de- 
creases as the temperature is increased. 

Steel spirals for reinforced concrete columns, SP 2535-6 
(Oct. 1, 1963), 10 cents. 

Establishes as a standard of practice in production, distri- 
bution, and use the sizes of steel spirals used for concrete 
column reinforcement in the building industry. The stand- 
ard establishes three standard sizes of steel bars and wire 
used for spirals, and includes a table listing the recommended 
pitches in multiples of 4s inch for spirals with core diameters 
from 11 to 39 inches inclusive for varying percentages of 
spiral reinforcement. Definitions of several terms applicable 
to steel spirals are also included. 

Vises (machinists’ and other bench-mounted vises), SP 2.229 
63, 10 cents. (Supersedes R229-48). 

Establishes a standard of simplified practice for the 
kinds, types, and sizes of bench-mounted vises. 
listed are currently in general use and 
considered to afford an adequate 
purposes. Accordingly, they are 
stock items. 


various 
The products 
demand, and are 
selection for ordinary 
recommended as standard 


Grading of a diamond powder in sub-sieve sizes, (S261-63, 
10 cents. (Supe rsedes in part CS123-49). 

Provides a nationally recognized standard for the sizing of 
diamond powder in sub-sieve (micron) sizes and a uniform 
method of determining conformance with the size designated. 
It serves as a common basis for understanding between 
purchasers and sellers as to the quality and particle size 
desired or supplied, and thereby should promote fair competi- 
tion. 


Absolute calibration of the National Bureau of Standards 
photoneutron source: III. Absorption in a heavy water 
solution of manganous sulfate, R. H. Noyce, Ek. R. Mosburg, 
Jr., S. B. Garfinkel, and R. S. Caswell, J. icl. Eng. 17, 
No. 7, 313-319 (1963). 

The neutron emission rate for the National Bureau of Stand- 
ards radium-beryllium (7, mn) standard source (NBS-I) 
has been redetermined by a relative comparison to an anti- 
mony-beryllium source which had been calibrated absolutely 
in a heavy water manganous sulphate solution through an 
indirect method involving 4%8= +7 coincidence counting of the 
induced Mn * activity. Correction was made for absorption 
of neutrons by other elements present in the bath. The 
emission rate of NBS-I was found to be 1-255 * 10° n/see 
as of June, 1961 with an uncertainty of 1 per cent. This 
measurement, when suitably averaged with previous deter- 
minations (DEJUREN, PapGerr and Curtiss, 1955; DeJUREN 
and CHIN, 1955) vields the best value 1-25; * 10° n/see with 
an uncertainty of 1 per cent. 


Determination of A ag oper pe X- ray photon flux and total 
beam energy, H. W. Koch and J. Pruit, Ed. L. C. Yugan 
and C. Wu, Nuel. Shae 5, ch. 2. 8. 3 508-553 (1963). 

The electron impact spectra of helium and ethylene are meas- 
ured at 50 volts primary energy and in the straight ahead 
direction. Under these conditions the spectra are analagous 
to the ultraviolet absorption spectra. It is concluded that 
previously reported deviations are due to the high angle col- 
lection geometry rather than the low voltage. 


Metallurgical microanalysis with the electron probe, J. R. 
Cuthill, L. L. Wyman, and H. Yakowitz, J. Metals AIME 
15, No. 10, 763-768 (Oct. 1963). 

The electron probe microanalyzer that was designed and built 
at the National Bureau of Standards is described, the func- 
tional design considerations involved are discussed, and the 
long-felt need for such an instrument in the metallurgical 
field is reviewed. 
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The electron probe microanalyzer makes possible for the first 
time the quantitative chemical analysis of microconstituents 
n situin an alloy, and provides a convenient method for deter- 
mining the concentration gradients alloy structure. 

The NBS electron probe microanalyzer employs a normal 
incidence, high-resolution optical system concentric with 
the axis of the electron beam It is equipped with two sean- 
ning X-ray spectrometers with a capability of analyzing for all 
of the e from Mg to U, if they can be suitably 
ined in a vacuum 

The application of the probe to two investiga- 
tions, representative of those carried out thus far, is discussed. 
(Additional features and modifications to the NBS electron 
probe, to be incorporated In the neal 1 
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Heating rate as a test of adiabatic calorimeters and the heat 
capacity of a alumina, I. 1). West, 7 Faraday Soc. 59, 
i, pt. 9, Sepl. 196 
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spectral irradiance, Stair, W. FE. 
Jackson tppl. Opt. 2, No. 11, 1151 
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Intercomparisons of the standard thermal-neutron flux density 
of the National Bureau of Standards, W. M. Murphey and 
J. Chi Proce S p \ on Detection, Dosimet and 
windardization, Int vional Atomic Enerqy Agence Llomu 
Res “stablishment, Harwell. England. Dec. 10 
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inte rpre tation of pH measurements in alcohol-water solvents, 

I . M. Paabo and R. A J. Phus. Chem. 

7, ] 196.5). 

pH numbers for solutions in 
iined from the e.m.f. of the 

electrode, salt bridge and 


tobinson, 


alcohol-water 
cell with hydrogen 


calomel reference elec- 


trode, standardized with aqueous buffer standards, are subject 
to no simple interpretation. Because of the indeterminate 
potential at the junction solution X (aleohol-H,O) |KCl 
bridge (aq.), they do not lie on a conventional scale of hydro- 
gen ion activity referred either to the aqueous standard state 
(pay) or to the standard state in the alcoholic medium 
(pay*). Values of # + log y Cl (where E; is the liquid- 
junction potential expressed in pH units and mc; is the pri- 
mary medium effect of chloride ion) were found to be as 
constant for different buffer solutions in alcoholic solvents of 
fixed composition as for the strictly aqueous medium. Inas- 
much as myc; at a given temperature is dependent only on 
solvent composition, the liquid-junction potential must there- 
fore nearly constant for a given solvent mdeium. 
Correction terms EF log m+H 6 have been calculated, and 
it is shown that pH—é6 closely approximates pay* under 
optimum conditions of measurement. An operational scale 
of pH* which is related to pay™ in the same way that pu is 
related to pay) is described. Values of pay* were obtained 
for 12 solutions in methanol-water solvents (0 to 68.1 wt. % 
MeOH) and for ethanol-water solvents (0 to 
100°7 EtOH 
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The error rate in a multiple-frequency shift system and the 
output signal noise ratio in a frequency modulation and a 
pulse-code-modulation frequency-shift system, H. Akima, 
Intern. Satellite London, Nov. 22-28, 
1963, pp. 305-309 Digest, London, 
England, 1963). 
Jused on the assumption of a fading-free 
tive white Gaussiun noise, the element 
an MFS (multiple-frequency-shift) system and the output 
SNR al-to-noise ratio) in KM frequency-modulation) 
and a PCM-FS (pulse-code-modulation-frequency-shift) 
tem are evaluated for wide ranges of system parameters. It 
is shown that the required intrinsic SNR for a given symbol 
rate in an MFS system can be reduced by increasing the 
number of frequencies in the keying. The possibility of 
improving the threshold of an FM system beyond that of a 
conventional one by frequeney modulating the carrier with 
sampled values and demodulation the modulated wave with 
a band-dividing demodulator is shown. A brief discussion 
on the threshold effects in the frequeney-lock and phase-lock 
FMI demodulators suggest that the threshold of these feed- 
back demodulators cannot be improved beyond that of a band- 
dividing one It is shown that the threshold in a PCM-FS 
svstem can be reduced by increasing the base in the coding. 
The comparison of a PCM-FS system with a band-dividing 
FM system shows the inferiority of the former to the latter, 
as the minimum-power criterion of system comparison 
is concerned, 
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A study of stress relaxation with finite strain, B. 
K. A. Kearsley, and L. J. Zapas, Trans. Soc 
391-410 (1963). 

Two simple types of equations appropriate to 
materials exhibiting elasticity are presented, one of a basic 
solid nature and one of a basic fluid nature. The predictions 
of the equations for a stress relaxation experiment are worked 
out and compared to the data from experiments on 
various elastomers. The fluid theory is shown to be most 
appropriate in a certain sense. 


Jernstein, 
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Experimental evaluation of a thallium beam ae uency stand- 
ard, R. KE. Beehler and D. J. Glaze, Pre 17th Annual 
Frequency Control Symp. 1963). 

\ thallium atomic beam frequeney standard has 
operation at the National Bureau of Standards since 
ber 1962 An evaluation of the system has been conducted 
in an effort to determine whether the use of thallium will 
result in a substantially improved atomic frequeney standard. 
The thallium system is a modified version of NBS I, the 
original Cs beam. 

The preliminary measurements show that, as expected, the 
existing uncertainties associated with the uniform Magnetic 
C field produce a corresponding uncertainty in the frequency 
measurements of only a few parts in 10*'%. Measurement 
precision has been determined from the reproducibility and 
statistical analysis of comparison data between the thallium 
standard and the USFS, NBS II. However, day-to-day 
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variations in the earlier comparisons are significantly larger 
than one would expect from the statistics. Possible sources 
of this difficulty along with improvements and the most 
recent data available are discussed. 

The best measured values for the frequency of the thallium 
resonance in terms of cesiumare given and compared with 
that published by Dr. Bonanomi’s group at Neuchatel, 
Switzerland. 


Some factors affecting the creep behavior of an age-harden- 
able alloy, W. D. Jenkins and W. A Willard, Trans. Quart. 
ASM 56, No. 3, 427-441 (Se pl. 1963). 

The effect of hardening, due to alloying and to precipitation 
preceding and during testing, on the creep behavior and the 
effect of creep on the hardening characteristics of a 94 percent 
nic kel—4.25 percent aluminum alloy (Duranickel) at 700, 900, 
and 1200° F are discussed. Alloying nickel with 4.25 percent 
aluminum increased the creep strength of the annealed metal 
by as much as 7-fold. This factor was further increased by 
aging prior to testing. However, the strength ratio (¢ 
Duranickel/s nickel) decreased as the temperature was de- 
creased or the creep rate increased. Although no single 
equation could describe the stress-creep rate or fracture time 
relations, the shape of these curves was dependent on both 
temperature and initial condition of the alloy. Furthermore, 
the second stage creep rates were inconsistent indices of the 
fracture times. The ductility of the alloy was signficantly 
affected by prior aging, test temperature and creep rate. 
Photomicrographs, illustrating types of fractures, intersecting, 
wavy and straight slip bands, twinning, recrystallized grains, 
and the density and distribution of (Ni;Al) precipitates, are 
presented to assist in the analysis of the modes of deformation 
occurring under differing test conditions. 


Measurement of weak magnetic fields by optical pumping 
methods, P. L. Bender, Bull. A mpere 9°, pt. I], 261-628 (1960) 
Three new types of magnetometer, the de alkali vapor magne- 
tometer, the self-oscillating alkali vapor magnetometer, and 
the helium magnetometer, all based on optical pumping, are 
compared. Practical limitations and performance in various 
environmental circumstances are considered. 


Use of transistors in Van de Graaff source leak controls, 
A. C. B. Richardson and L. Costrell, Nuclear Instr. Methods 
Letter to Editor 24, No. 1, 129-130 (July 1963). 

The feasibility of using transistors in the high voltage ter- 
minals of pressurized accelerators has been demonstrated in 
a circuit developed for adjustment of the heater current of 
palladium leaks for the ion source for Van de Graaff accelera- 
tors. Improved resolution and stability, together with 
trouble free operation have been achieved. 


Recent 
R. S. Caswell, R. H. 
Mosburg, Jr., (Proc. Symp. Neutron 
and Standardization, International 
Atomic Energy Research Establishment, Harwell, England, 
Dec. 10-14, 1963), Book, Neutron II, 547-556 
(International Atomic Enerqy 1963). 

A new absolute measurement of the neutron emission rate 
of the U.S. national standard neutron source, a Ra-Pe(y,n) 
source has been using a manganous sulfate bath primarily 
of heavy water, which is less dependent on values of neutron 
absorption cross sections than previous bath methods. To 
avoid production of neutrons in the heavy water from D(y,n) 
reactions, an Sb-Be(y,n) source was absolutely calibrated 
in the heavy water bath, and in turn the standard source was 
calibrated by a relative comparison to the Sb-Be(y,n) source 
in a light water bath. The result of this measurement 
was a neutron emission rate of 1.24, 10° neutrons/sec, 
which when combined with previous measurements gives a 
best value of 1.25. * 10° neutrons/see. The estimated 
uncertainty is 1 per cent. 

Measurements of the growth of neutron emission rate in 
Pu-Be(a, n) sources have shown these sources to have varying 
rates of growth of neutron emission up to about 2 percent 
year. For use as standards, information on isotopic composi- 
tion or on neutron emission rate growth is required. Am ?4!- 
Be(a,n) sources offer promise of providing a standard source 
with characteristics similar to the Pu-Re sources, but without 
the growth uncertainties 


developments in neutron source standardization, 


Noyee, 8. B. Garfinkel, and E. R. 
Detection, Dosimetry 


Atomic Energy Agency, 


Dosimetry 


Agency, Vienna, 


Some physical properties of organic denture base materials. 
J. B. Woelfel, G. C. Paffenbarger, and W. T. Sweeney, /. Am, 
Dental Assoc. 67, 489-504 (Oct. 1963). 

Pertinent physical properties were measured of 12 different 
types of organic denture base materials used in the fabrication 
of technic and clinical dentures which had been under observa- 
tion for from three to four years. The cold-curing acrylic 
resins, the weakest and the most flexible did not fracture or 
distort in use. The polystyrene resin, the strongest and 
least flexible fractured in two out of five upper dentures. 
It is surmised that these fatigue-like fractures were caused by 
a high residual stress imparted during molding. The repaired 
resins had 40-71% of their original strength. The Izod 
impact strength ranged from 0.17 (heat-curing acrylic resins) 
to 1.65 ft.-lb. per inch of notch (glass fiber filled acrylic resins). 
Indentation tests showed the polystyrene to be the hardest. 
The coefficients of linear thermal expansion of ‘“wet’’ speci- 
mens run at 20—37° C. ranged from 42.7 10-*/°C for glass- 
fiber-filled acrylic resins to 98.3>/ 10-*/°C for a heat-curing 
acrylic resin copolymer of the powder-liquid type. Both 
laboratory and clinical tests are needed in the evaluation of 
denture base resins. 


Dependence of power output of a gas laser on the length and 
rate of excitation of the discharge, J. A. White, Appl. Phys. 
Letters 3, No. 7, 107-109 (Oct. 1, 1963). 

The gain and emittance of a gas laser are calculated from 
first principles for the case that the atoms in the gas interact 
appreciably with only a single Doppler-shifted component of 
the radiation. Good agreement is obtained with the experi- 
mentally observed output of a laser operating simultaneously 
in several modes. Although the theory does not take into 
account interference effects when the same atoms are per- 
turbed appreciably by different modes of oscillation or dif- 
ferent Doppler components of the same mode, it is otherwise 
applicable for arbitrarily high field intensities. 


The adsorption of methane and nitrogen on silica gel, syn- 
thetic zeolite and charcoal, A. J. Kidnay and M. J. Hiza, 
J. Phys. Chem. 69, 1725-1727 (1963). 

The data obtained in this study provide a characterization of 
the temperature dependence of methane adsorption near the 
boiling point of liquid nitrogen and a comparison of the 
adsorption isotherms of methane and nitrogen on three 
commercial adsorbents. Isosteres of nitrogen on silica gel 
and methane on silica gel, charcoal and a synthetic zeolite 
from 70° to 82°K are presented; in addition isotherms for 
methane and nitrogen on all three adsorgents at 76°K are 
presented. The were used to calculate 
thermodynamic properties of these systems. 


isosteres some 


Galvanostalametry, a new technique based on the negative 
pressure of liquids, for investigating electrechemical phe- 
nomena at an electrode, J. Sligh and A. Brenner, /. Flec- 
trochem. Soc. 110, No. 11, 1136-1142 (Nov. 1963). 
A procedure based on the negative pressure of water has 
been used as a sensitive indicator for the study of electro- 
chemical phenomena. The apparatus consists of a vertical 
glass tube, closed at the upper end, which is evacuated and 
filled with an electrolyte. The latter remains suspended in 
the column in a metastable state of tension. The column of 
electrolyte was dropped by producing a minute amount of 
electrolysis between 2 electrodes sealed in the top of the 
column. The time or current required to drop the column 
(abbr. t.d.c.) was used for making the following experiments 
and/or observations: 
A. The decomposition potential of water was determined 
by observing the lowest applied voltage required 
t.d.e. 

The time required t.d.c. (on passage of a constant 
current) was proportional to the concentration of 
iodide ion in the electrolyte. 

The current required t.d.c. was proportional to the 
area of the indicator electrode, regardless of the 
irregularity of its shape. 

The formation of the gas required t.d.c. occurred in 
less than 15 microseconds after closing the circuit. 
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Electrical resistance-strain characteristics of thin evaporated 
metal films, R. L. Parker and A. Krinsky, Appl. Phys. 
34, No. 9, 2700-2708 (Se pt. 19638 

The electrical resistance-strain coefficients of thin evaporated 
films of aluminum, gold, cobalt, nickel, palladium, platinum, 
antimony, and tellurium have been measured. The coeffi- 
cients strain-sensitivity coefficient or “gauge-factor” 4 
depend on the specific resistance R, of the films. For 
relatively thick films (small R tends to approach the 
bulk value for relatively thin films (large R,), y may be 
many times greater than for films of intermediate thick- 
minimum value The resistance change was 
found to be proportional to strain for all the above films except 
those of Ni and Sb. possible mechanisms for the 
observed behavior are it is suggested that free-path 
effeets are predominant in films of intermediate thickness, 
and that the tunneling conduction mechanism plays an 
important role in producing the high strain sensitivity of very 
thin (high PR films of separated particles. 


hess has a 


Some 


discussed 


A 100 megapulse per second binary counter with impedance 
soaring, Ss. B. Geller and P. A. Mantek, JEEE T, 
Electronic Computers EC-12, No. 5, 568 (Oct. 1963 

This Letter 100 meg apulae-per-aceons triggerable 
flip-flop or binary counter circuit. It employs transistor- 
tunnel unique trigger steering 
impedance states in the 


rans. 
aeseripes a 


diode circuitry 


mechanism that depe nds 


and has a 
upon the 


circult 


Line shapes, saturation behavior, and temperature studies in 
the nuclear resonance of nickel, R. L. Streever and L. H. 
Bennett, Phys. Rev. 1381, No. 5 (Sept. 1, 1963 
The 1 magnetic Ni®! has been studied in 
natural nickel (1.25°% abun- 
( i ind hed nickel (99.90) abundant in 
Ni were used The signal with a — 
tive to the a combination of a bsorption 
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iration behavior an estimate of 
applied rf level for nuclei in domain 
made An additional broadening 
observed in the enriched nickel, particularly at 
Application of Redfield’s theory of saturation at high 
leads to the hat the additional broadening 
direct exchange between nuclei. 
isurements from room temperature at 4.2°K 
i marginal both the enriched and 
ckel are \ temps rature-depende nt 
ienev shift the enriched and unenriched 
nickel which appears to be a maximum at about 200°K where 
the enriched nickel frequency is about 30 ke/see higher. From 
temperature to 536°K, frequency measurements made 
super-regenerative circuits are reported. A comparison 
the temperature dependence of the resonance frequency 
the saturation magnetization is made. 
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Impedance of Leclanché cells and batteries, R. J. 
H. J. DeWane, Electrochem. Soc. 110, No. 11, 
Vov. 196 


Bro dd 
1091 


and 
109) 


The impedance, 


performance on 


open circuit voltage, 
standard tests, an 
letermined for a number of diff 
iced by various manufactu 
impedance in the 


short circuit current, 
d the pulse resistance were 
ferent size Leclanché cells 
rers. The interpretation of 
light of modern perturbation theory 
ndicated the existence of one process at the negative zine 
lectrode and two processes at the positive manganese 
lioxide electrode for all Leclanché cells assuming ideal behav- 
electrode process The magnitude of the con- 
tributions to the total impedance and the relaxation time for 
each of the electrode processes in the cell were determined. 
The effect of varving dioxide to carbon black 
weight ratios on the impedance of laboratory-prepared cells 
determined Also, the impedance of fresh cells and 
batteries was determined at several stages during discharges 


1O for each 


manganese 


was 


on various standard tests, and, in some cases, after shelf stands 
of up to two years. The extrapolated value of the impedance 
at infinite frequency, Ro, was used to estimate the residual 
capacity of cells after shelf stands of one and two years. 
Reo was found to be identical with the internal resistance of 
the cell determined by pulse techniques. 


Absolute scale of oscillator strengths, ( 
Corliss, Monthly Notices Roy. 
(1963). 

Absolute f-values assembled at the University of London 
Observatory (U.L.O.) using the f-sum rule, and at the 
National Bureau of Standards (N.B.S.) using experimental 
values are compared. The systematic logarithmic difference 
is U.L.O.—N.B.S. 0.19 0.04 dex. From the discussion 
it is concluded that either method may be used for obtaining 
extensive absolute f-values with an accuracy of about + 0.2 
dex. 


Allen and C, H. 
126, No. 1, 37 40 


Astron. Soc. 


Relation of partial (110) pole figures to thickness and micro- 
structure of electrodeposited copper, F. Ogburn and C. J. 
Newton, + Electrochem. Soc. 110, No. 11, 1148-1150 (Nov. 
1963). 

The distributions of (110) poles near the surfaces of two copper 
electrodeposits are given. The concentration of (110) poles is 
maximum for y=0° and minimum for y=30° where y is the 
angle between the (110) pole and the normal to the specimen 
— The specimen that was 1700u thick showed a much 
greater concentration near y=0° than the thinner deposit of 
l * . A cross-section micrograph of the thicker deposit shows 
many twinned planes near perpendicular to the deposit 
surface. 


Regulation * are gauge emission current to better 
than 0.05°; Yee and R. J. Carpenter, Rev. Sci. Instr. 
34, No. 10, vies ye (Oct. 1963). 
A control svstem is described for 
limited emission current from the filament of an ionization 
gauge tube. A d-e transistor circuit is used with a lamp and 
a photosensitive resistor to provide electrical isolation be- 
tween the control and the controlled circuits. Emission cur- 
rents from 10 mA to 0.05 mA can be maintained with a sta- 
bility of better than 0.05°% with the system. 


regulating the temperature 


Errors in “‘drop”’ calorimetry due to sample container transi- 
tions, D. C. Ginnings, J. Phys. Chem. 67, 1917-1918 (1968). 
An examination was made of three sets of published enthalpy 
data obtained by using the “drop”? method and container 
materials having transitions It is concluded that in all 
three cases, small but significant errors resulted from the use 
of the container materials with transitions. 


Energy dependence of proportional-counter fast-neutron do- 
simeters, R. 8. Caswell, W. B. Beverly, and V. Speigel, Jr., 
(Proc. S Neutron Detection. Dosimetry and Standardiza- 
fron, tomic Enerqu, Atomic Enerqy Research Establishment, 
Harwell, Enaland, Dec. 10-14. 1962), Book, Neutron Dosime- 
try I, 227-237 (July 1963 

Proportional counter dosimeters with polyethylene walls and 
ethylene gas filling are widely used for 
ments in the presence of gamma 
differentiated by pulse-height 
larger than a given bias are accepted and pulse-height inte- 
grated to determine the neutron dose. In the present work, 
the energy response of a fast neutron dosimeter with a known 
bias was measured using known fluxes of monoenergetic 
neutrons between 0.1 and 4 Mev. The pulse-height distri- 
butions were also observed. An approximate calculation has 
been made which predicts reasonably well the energy depend- 
ence of this type of dosimeter. More detailed calculations 
should be able to predict quite accurately the response of 
this type of detector. 


jmp. 


neutron dose measure- 
rays, the radiations being 
discrimination. Only pulses 


Projection X-ray microscopy 
Newman and §. H. Greenfeld, 
893-898 (Nov. 1963). 

Point projection microradiography has been used to elucidate 
the relationship of components in structures containing as- 
phalt. Soft X-rays, in the 8 A range, produced good contrast 


of roofing materials, S. B. 
Mater. Res. Std. 3, No. 11, 
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in a number of specimens of roofing systems. Low magnifi- 
cations were adequate for indicating defects and advantages of 
the materials for the intended applications. The method 
provides a unique tool for obtaining such information. 


A redetermination of some optical properties of calcium 
fluoride, I. H. Malitson, Appl. Opt. 2, No. 11, 1103-1107 (Nov. 
1963) 

This paper gives refractive-index values at 24°C for a syn- 
thetic and a natural prism of calcium flouride. Measure- 
ments were made by the minimum-deviation method at 46 
calibrated wavelenghts from 0.23 yw in the ultraviolet to 9.7 pu 
in the infrared. The indices of the synthetic prism were 
fitted to a three-term Sellmeier dispersion formula of the form: 
n?—1 = DAjd?2/(A\27—A?? Dispersive quantities were com- 
puted which estimate the expected relative dispersion, 
chromatic aberration, and resolution of CaF, as a function 
of wavelength. Valeus of dn/dt are also given. The results 
of this work are compared with previously reported data. 


Calorimetric method for measuring polyester degradation 
due to weathering, V. f. Gray and J. R. Wright, J. Appl. 
Polymer Sci. 7, 2161-2174 (1968). 

Chemical changes occurring in thermoset polyesters during 
exposure to solar and carbon-are radiant energy were deter- 
mined quantitatively by a color reaction with a diamine, 
N,N-dimethyl-p-phenylenediamine, in a methanol-benzene 
solution. The amount of diamine reacting with a specimen 
and the color change of the treated specimen, as measured by 
differential reflectance, varied directly with exposure time. 
Three commercial polyesters (fiberglass reinforced) exposed 
under identical conditions had widely different degradation 
rates as measured by the amount of N, N-dimethyl-p-phenyl- 
enediamine used per specimen The relative stability of the 
three materials was independent of irradiation source. The 
color of the treated specimen was dependent on exposure 
source, 


Optical T-bench method for measurement of optical path 
difference, F. k. Washer and W. R. Darling, SPJE J. 1, 
211-217 (Sept. 1963 


A visual optical T-bench method for measurement of optical 


path difference n, in units of wavelength is described. The 
method employs an optical T-bench equipped with nodal slide, 
a movable double-slit system, and an angle measuring tele- 
Determinations of longitudinal spherical aberration 
in millimeters and optical path difference in units of wave- 
length are upon measurements of angular deviations 
in small selected regions of the collimated beam emergent 
from the lens under test. The theory of the method is pre- 
sented together with a description of the apparatus used 
and the technique of measurements. tesults of measure- 
ment on a single lens are included. 


scope 


based 


Programmed maneuver-spectrum fatigue tests of aircraft 
beam specimens, L. Mordfin and N. Halsey, ASTM STP 
Symp. Fatigue Tests of Aircraft Structures, No. 338, 251-272 
(Sept. 1963). 

Constant- and variable-amplitude fatigue tests were con- 
ducted on 7075-T6 aluminum alloy built-up beam specimens 
in bending to study certain aspects of spectrum fatigue. 
Existing cumulative fatigue damage theories were found in- 
adequate for the dual purpose of predicting spectrum fatigue 
life and evaluating the effects of individual load levels in a 
spectrum. However, the observed spectrum fatigue prop- 
erties were consistent with the measured effects of preloading, 
overloading, and underloading. 

Preloading was beneficial when the preload was considerably 
greater than the subsequent fatigue loads, and when the pre- 
load was applied in the same direction as the fatigue loads. 
When the preload and the fatigue loads had opposite signs, 
a reduction in fatigue life resulted. These effects were magni- 
fied by the periodic application of overloads. 

The periodic application of repeated low loads was ineffectual 
when applied in the same direction as the fatigue loads. 
When the spectrum included low loads of negative sign, serious 
reductions in life resulted. 

The most valuable finding, perhaps, was that the fatigue be- 
havior of the beam specimens was similar to that of certain 


full-scale aircraft structures. It was also found that the 
fatigue properties of the beam specimens were rather insensi- 
tive to variations in testing technique and material. 

A pneumatic testing machine and an automatic programmer 
that were developed for this investigation are described. 


Differential dielectric apparatus for determining water added 
to solvents, W. C. Wolfe, Anal. Chem. 35, 1884-1887 (Nov. 
1963). 

A simple and rapid technique 
small changes in water 
(2-ethoxyethyl) ether 


is described for determining 
concentration in p-dioxane or bis 
(diethylene glycol diethyl ether). The 
moisture content of solids can be determined by extraction, 
followed by this procedure. Dual matched dielectric measur- 
ing cells are filled with liquids to be compared and capacitance 
differences are measured with a commercial Q-meter. The 
cells are calibrated with the solvent to be used and solutions 
containing known concentrations of water in the same solvent. 
Values for change in dielectric constant with change in water 
concentration (weight percent), de/dC, are calculated from 
experimental data and from literature values. The tempera- 
ture coefficient of de/dC is sufficiently = ill so that precise 
temperature control is unnecessary. Calculations of dipole 
moment from the data support Sie observations that 
water behaves in a non-associated manner in dilute 
in p-dioxane. 
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An interpretation of rapid changes in the phase of horizontally 
polarized VLF waves recorded at night over a short path 
in the southwestern United States. E. E. Gossard and 
M. R. Paulson. 

Precise phase and amplitude measurements on VLF ign als 
propagated through the Arctic zone. F. H. Reder, C. J. 
Abom, and G. M. R. Winkler. (See above Shamans.) 
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On the long term phase stability of the 19.8 ke/s signal trans- 
mitted from Hawaii, and received at Boulder, 
A. H. Brady. (See above abstracts.) 

Oblique propagation of 
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transmission of the atmosphere to solar radiation, 


D. M. Gates and 
(Sept. 1963) 

Synoptic variations and vertical profiles of large-scale iono- 
spheric irregularities, R. G. Merrill, 
N. J. Roper, J. Geophys. Res. 68, 
(Oct. 1, 1963). 

Excitation of the red and green coronal lines, C. Pecker and 
R. N. Thomas, Ann. Astrophys. 25, No. 2, 100-108 (1962). 

Some topics in quantum statistics. The Wigner function and 
transport theory, H. Mori, I. Oppenheim, and J. Ross, 
Book, Studies in Statistical Mechanies, ed. DeBoer and 
Uhlenbeck, I, pt. C, 218-298 (North Holland Publ. Co., 
Amsterdam, The Netherlands, 1962). 

Normal congruence subgroups of the (> 
Newman, Bull. Am. Math. 
No. 5, 719-720 (Sept. 1963 

Thermodynamic study of the thorium phosphide with a mass 
spectrometer, J. Efimenko and K. A. Gingerich, 
Symp. Thermodynamics of Nuclear Materials, 
May 21-25, 1962, pp. 477-486 (International 
Energy Agency, Vienna, Austria, Sept. 1962). 

Radiation impedance of a source near reflectors, R. V. Water- 
house, J. Acoust. Soc. Am. 35, No. 8, 1144-1151 (Aug. 1963). 

Low energy levels of neutral cerium (Ce 1), W. C. Martin, J. 
Opt. Soc. Am. 53, No. 9, 1047-1050 (Sept. 1963). 

The energy environment in which we live, D. M. Gates, Am. 
Scientist 51, No. 3, 327-348 (Sept. 1963 

Evidence for field-aligned ionization irregularities between 
100 to 100 km above the earth’s surface, T. E. VanZandt, 
W. Calvert, R. W. Knecht and G. B. Goe, Third Intern. 
Space Sci. Syvmp. Proc., Washington, D.C., May 2-8, 1962, 
pp. 271-273 (John Wiley & Sons Inc., New York, N.Y., 
1962). 

Studies in the geomagnetic micropulsation frequency range, 
W. H. Campbell, IEEE Intern. Conv. Record 11, pt. 9, 
25-28 (Mar. 25-28, 1963 

Determination of the electronic energy levels of molecules by 
low energy electron impact sp ctroscopy, J A. 
and 8S. R. Mielezarek, J. Phys. 
1963). 
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